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FAR  INFRARED  ABSORPTION  BY  SMALL  PARTICLES 

By 

Young  Hoon  Kim 
May  1986 

Chairman:  David  B.  Tanner 

Major  Department:  Physics 

The  properties  of  small  metal  particles  of  the  order  of  100  A 
size  have  been  studied.  In  particular,  the  anomalously  large  far 
infrared  absorption  by  small  particles  has  been  a puzzling  problem 
since  its  discovery  in  1972.  Even  though  many  experimental  and 
theoretical  efforts  have  been  made  to  understand  this  peculiar 
phenomenon,  no  clear  answer  has  been  made  yet.  In  order  to  understand 
the  underlying  physics  of  the  anomalous  far  infrared  response  of  a 
small  "metallic"  system  whose  geometrical  dimension  is  smaller  than 
any  characteristic  length  scale,  a question  about  the  fundamental 
nature  of  a small  particle  was  raised  and  investigated  both 
experimentally  and  theoretical ly  in  this  work. 

For  the  measurements,  small  metal  particles  (A1 , Ag,  or  Au) 
were  mixed  with  insulating  particles  (KCl  or  paraffin),  with  the  small 
particle-insulator  mixture  being  then  compressed  into  a pellet. 
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Transmission  electron  microscopy  of  ion-milled  or  chemically  polished 
samples  showed  the  random  nature  of  the  small  particle  composite. 

The  far  infrared  through  visible  frequency  optical  properties  of 
samples  having  metal  fractions  between  0.0001  and  0.032  were  measured. 
A persistent  anomalously  large  far  infrared  absorption  was  observed, 
which  could  not  be  explained  by  existing  theories. 

In  the  mid-infrared  frequency  range,  the  absorption  was 
smaller  and  could  be  fit  to  classical  theories  by  assuming  an 
unexpectedly  short  electronic  mean  free  path.  An  extremely  broad 
surface  plasmon  resonance  was  observed  for  Ag-KCl  and  Au-KCl  composite 
system.  Additionally,  different  energy  band  structure  from  the  bulk 
state  was  found  for  aluminum  small  particles. 

Theoretical  work  included  calculations  of  the  dielectric 
function  of  a small  particle  at  frequencies  near  the  size-quantization 
energy  gap.  These  calculations  showed  that  the  particle  behaved  as  an 
insulator  rather  than  as  a small  solid  metal. 
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CHAPTER  I 
INTRODUCTION 

This  dissertation  explores  the  far-infrared  (FIR),  the 
mid-infrared  (MIR),  and  the  optical  properties  of  small  "metallic" 
particles  experimentally  and  theoretically.  In  this  work,  three 
different  kinds  of  small  particles  (A1 , Ag,  and  Au)  were  studied.  These 
small  particles  were  dispersed  in  an  insulating  medium  (KCl  or  paraffin) 
for  experimental  measurements  of  the  electromagnetic  response. 

One  might  expect  that  the  properties  of  a small  particle  would 
be  different  from  those  of  an  infinite  metallic  system  when  the 
geometrical  size  of  a metal  is  smaller  than  the  physical  characteristic 
length  scale.  Also,  one  might  think  some  unusual  physical  properties 
arise  from  the  smallness  of  the  size  of  the  particles.  This  aspect  will 
be  extensively  studied  and  discussed  in  this  work. 

Statement  of  the  Problem 

Metals 

A metal  is  a type  of  condensed  matter  with  high  electrical 

conductivity.  A great  detail  of  metal  physics,  both  the  static  and  the 

dynamic  properties  of  the  electrons  such  as  heat  capacity,  magnetic 

susceptibility,  electron  spin  resonance,  and  electromagnetic  properties, 
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etc.,  is  now  understood.  ’ 

When  isolated  metal  atoms  are  brought  together  to  form  a solid, 
the  valence  electrons  interact  with  their  counterparts  on  neighboring 
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atoms.  As  a result  of  this  interaction  the  electrons  are  no  longer 
effectively  associated  with  any  particular  constituent  atom.  In 
contrast,  the  core  electrons  are  generally  unaffected  by  the 
condensation  and  remain  bound  to  individual  nuclei.  Consequently,  we  can 
view  a solid  crystalline  metal  as  having  two  structural  units.  One  is  the 
group  of  positive  ions  which  are  on  the  regular  periodic  lattice  sites 
and  which  usually  are  treated  as  a positive  jellium  background  in  a 
theoretical  calculation.  The  other  is  the  group  of  conduction  electrons, 
which  are  distributed  throughout  the  interionic  space  and  which  often  are 
modeled  as  an  electron  gas  or  a Fermi  liquid.  At  zero  temperature  these 
conduction  electrons  fill  all  available  energy  levels  up  to  a maximum 
energy  which  is  called  the  Fermi  energy.  The  energy  distribution 
function  of  the  electrons  has  a sharp  discontinuity  at  the  Fermi  energy 
regardless  of  any  interactions,  although  interactions  such  as  impurity 
scattering  causes  a smearing  of  the  discontinuous  step  in  their  momentum 
distribution  function.  Therefore,  we  can  characterize  a metal  as  a 
system  which  has  a Fermi  surface  and  easily  accessible,  vacant  excited 
energy  states  available  to  the  conduction  electrons. 

The  existence  of  the  conduction  electrons  leads  to  the 
interesting  physical  properties  of  a metal  such  as  the  conduction  electron 
energy  band  which  consists  of  many  energy  levels  and  a collective 
excitation  (plasma  oscillation)  of  the  free  electrons.  The  origin  of 
the  discrete  energy  levels  in  which  the  electrons  can  be  excited  in 
response  to  the  external  influences,  such  as  external  electromagnetic 
fields  or  temperature  gradient,  lies  in  the  quantization  of  the  kinetic 
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energy  of  electrons.  Since  there  are  about  10  electrons/cm  in  a 
typical  metal,  their  average  energy  level  spacing  would  be  the  order  of 
10'^*^  eV,  which  essentially  represents  the  continuum  states. 

In  order  to  study  the  physics  of  the  conduction  electrons,  an 
important  quantity  to  understand  is  the  electrical  conductivity  (or 
dielectric  function)  of  a metal.  If  there  was  no  lattice  excitation 
(phonon)  and  no  lattice  defects  in  an  entirely  pure  metal,  there  would 
be  no  scattering  of  the  electrons  by  the  structure  and  the  conductivity 
a (oj)  of  the  system  would  be 

0 (o))  = I 6 (oi)  0-1) 

where  n is  the  electron  number  density,  m is  the  electron  mass,  and  w is 
the  frequency  of  the  applied  fields.  Equation  (1.1)  describes  a system 
in  which  the  lifetime  of  the  electrons  is  infinity.  However,  because  in 
actual  metals  the  conduction  electrons  interact  with  each  other,  with 
the  ions,  and  especially  with  the  impurity  atoms,  the  lifetime  of  a free 
electron  between  collisions  becomes  rather  short,  leading  to  a finite  dc 
conductivity  and  the  optical  conductivity  as 


a (u) 


1 + iwT 


(1.2) 


Here  a = and  x is  the  relaxation  time,  which  is  assumed  to  be 
0 m 

constant. 

However,  one  difficulty  in  the  free  electron  gas  approximation  is 
that  a free  electron  gas  cannot  absorb  a photon  because  an  electron 
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cannot  conserve  both  momentum  and  energy  by  absorbing  the  photon  alone; 
therefore,  there  is  no  mechanism  for  the  loss  of  electromagnetic  energy 
for  a metal.  Instead,  a metal  actually  absorbs  some  small  amount  of 
radiation  energy  in  the  FIR  region  where  the  electrons  in  metals  show  an 

3 

outstanding  free  electron  gas  nature.  Holstein  showed  that  there  exists 
two  dissipative  processes  in  a metal:  One  is  due  to  the  presence  of  a 

surface  which  enables  the  electrons  to  acquire  an  energy  from  the 
oscillating  external  electric  field  when  they  pass  through  the  skin 
depth  layer,  and  then  the  absorbed  energy  is  converted  into  heat  through 
the  diffuse  collision  with  the  surface.  The  other  is  a bulk  process 
which  is  known  as  the  Holstein  process.  The  conduction  electron  absorbs 
an  incident  photon  and  emits  a phonon  simultaneously.  Consequently,  a 
metal  does  absorb  the  FIR  photon  because  the  collision  with  the  surface 
or  the  Holstein  process  allows  the  conservation  of  both  energy  and 
momentum  for  an  electron.  But,  these  realistic  contributions  modify  the 
relaxation  time  of  an  electron  to  be  weakly  frequency  dependent,  modify 
the  Drude  behavior  to  some  extent. 

Now,  what  will  happen  if  we  reduce  the  geometrical  dimension  of 
a metal  to  a scale  which  is  smaller  than  any  physical  characteristic 
length  such  as  the  electron  mean  free  path  in  a metal,  penetration 
depth  of  fields,  the  wave  length  of  the  external  fields,  and  the  electron 
coherence  length,  but  still  large  enough  to  be  called  a small  solid 
metal ? 
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Small  Particles 

Because  the  size  of  the  metal  is  smaller  than  the  electron  mean 

free  path,  the  conduction  electrons  now  can  see  the  surface.  In  other 

words,  the  translational  symmetry  in  the  motion  of  electrons  is  broken. 

One  of  the  consequences  of  the  broken  symmetry  is  that  the  surface  plays 

a direct  quantum  mechanical  role  and  actually  determines  the  eigenstates 

of  electrons  bound  within  the  finite  volume  of  a particle.  If  we  model 

a small  particle  as  a spherical  box  filled  with  electron  gas  and  assume 

that  the  smooth  and  regular  surface  serves  as  an  infinite  potential 

barrier,  the  corresponding  eigenstates  would  be  highly  degenerate,  and 

1/3 

the  average  energy  gap  AE  would  be  the  order  of  Ep/N^  where  Ep  is  the 
Fermi  energy  and  is  the  number  of  electrons  in  the  system.  But,  as 
proposed  first  by  Kubo^  the  random  irregularity  of  the  surface  is 
sufficient  to  break  the  degeneracy  and  to  randomize  the  distribution  of 
levels  (a  simple  illustration  is  presented  in  Appendix  A).  Therefore, 
the  corresponding  energy  gap  would  be  the  order  of  Ep/N^  if  all  the 
degeneracy  is  lifted  and  the  distribution  would  follow  a Poisson 
distribution.  As  shown  in  Figure  1.1,  a model  of  a small  particle  with 
an  FCC  structure  reconstructed  from  the  Field  Ion  Microscope  analysis  of 

5 

tungsten  is  indeed  rough  to  atomic  scale  and  indirectly  justifies  the 
choice  of  the  random  distribution  of  the  energy  levels.  However,  the 
correlation  between  the  energy  levels  in  such  a random  system  imposed 
an  importance  of  the  statistics  of  new  classes  of  ensemble  which  was 
hypothesized  by  Dyson^  and  was  first  applied  to  the  small  particle 
problem  by  Gor’kov  and  Eliashberg^  (see  Appendix  B). 


Figure  1.1  Small  particle  model 

(a)  A model  of  a small  particle  reconstructed  from  the  Field 
Ion  Microscope  analysis.  (Reference  5) 

(b)  A Field  Ion  Micrograph  of  Tungsten  tip  (Reference  5); 
the  center  axis  is  [Oil] 
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Another  consequence  of  the  smallness  of  a metal  particle  is  the 

broken  synmetry  in  the  Bloch  state  at  the  surface  which  has  been  ignored 

in  the  experiment  and  many  body  theories  for  a bulk  metal  because  of  the 

negligible  surface  to  volume  ratio.  As  the  size  of  the  system  gets 

smaller,  however,  this  ratio  no  longer  is  negligible.  Because  the  ratio 

of  the  number  of  the  surface  atoms,  N^,  to  the  number  of  the  total  atoms, 

-1  /3 

N,  for  a small  particle  is  proportional  to  N , the  role  of  the  surface 

imperfections  (vacancies  and  impurities)  which  localizes  electrons  at  the 

surface  and  forms  of  the  localized  band  between  the  conduction  and  the 

8 9 

valence  bands  of  a metal  ’ should  be  appreciable  experimentally  for  a 
small  particle  as  a result  of  the  quantum  size  effect. 

Lower  Bound  to  the  Smallness 

It  is  instructive  to  consider  a lower  limit  to  the  size  of  a 
small  particle  to  retain  the  metallic  characteristic  (free  conduction 
electrons).  Obviously,  as  the  size  of  a particle  decreases,  the 
metallic  nature  of  a particle  will  disappear  and  the  particle  will 
eventually  recover  the  molecular  behavior.  In  other  words,  we  can  set 
a criterion  to  determine  the  metal -insulator  (molecule)  transition  of 
a small  particle  by  considering  two  energy  scales  in  a similar  way  to  that 
done  by  Thouless^^:  One  is  the  gap  energy  (AE)  due  to  the  level 
quantization 

" 3Ng 


(1.3) 
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and  the  other  is  the  energy  (6E)  associated  with  the  diffusion  time  tq 
6E  = — (1.4) 

Now,  when  the  gap  energy  wins  over  6E  (i.e.,  AE/6E  >1),  the 
electrons  are  no  longer  free  electrons.  Also,  the  electrons  are  in 
the  extended  states  when  AE/6E  < 1.  Therefore,  at  the  transition  point 
(AE76E  = 1),  we  have  a critical  particle  radius  a^  as 


'Mo 

Zimh 


1/3 


(1.5) 


where  n is  the  electron  density.  For  example,  a(.(Al)  = 5 ^ with 
Tp  = 0.8  X 10”^^  sec,  n = 18.06  x 10^^,  and  Ep  = 11.63  eV  for  aluminum 
and  for  silver,  a^(Ag)  = 9.7  A with  = 4.0  x 10  sec,  n = 

90 

5.85  X 10  , and  Ep  = 5.48  eV.  Here  the  for  the  bulk  state  has  been 

used  because  we  reduce  the  geometrical  dimension  of  a particle  from  its 
bulk  state. 

Also,  this  criterion  can  be  compared  with  the  loffe-Regel 
metal-insulator  transition  criterion^^ 


kpa  < 2tt  (1.6) 

which  gives  a(Al)  = 3.59  A and  a(Ag)  = 5.24  A. 

Therefore,  the  argument  given  in  this  section  assures  that  we  are 
dealing  with  small  "metal"  particles  and  justifies  the  theoretical  model 
of  electron  gas  in  a restricted  geometry. 
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Anomalous  Far-Infrared  (FIR)  Absorption 

It  has  been  understood  that  the  classical  theory  of  scattering 
and  abosrption  of  electromagnetic  wave  by  a small  particle,  which  is 
known  as  the  Mie  theory,  is  almost  complete.  One  of  the  most  successful 
achievements  of  the  classical  scattering  theory  is  the  prediction  of  the 
surface  plasmon  mode  in  a small  particle  which  originates  from  the 
limitation  of  the  long-range  correlation  of  the  electrons  by  the 
surface  (i.e.,  a collective  oscillation  caused  by  a restoring  force 
from  the  induced  charge  on  the  surface).  This  is  also  known  as  the  Mie 
or  Maxwel 1 -Garnett  resonance. 

From  the  classical  point  of  view,  introducing  the  surface  means 
that  the  scattering  of  electrons  by  the  surface  should  be  included.  This 
can  be  done  by  adding  a surface  scattering  contribution  to  the 
relaxation  time,  t,  of  the  electrons  in  the  Drude  formula  given  by 
Equation  (1.2) 


where  t.  is  the  relaxation  time  for  the  electrons  in  a bulk  metal, 
b 

= a/Vp,  a is  the  radius  of  the  particle  and  Vp  is  the  Fermi  velocity. 
Because  the  Drude  dielectric  function  with  the  modified  relaxation 
time  explains  qualitatively  the  observed  small  particle  effect  in  the 

IP 

visible  frequency  range,  it  had  been  believed  that  the  same  is  true 
in  the  FIR  range  as  shown  in  Figure  1.2  until  the  anomalous  FIR 
absorption  by  small  metallic  particles  which  shows  a contradiction  to  the 
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Figure  1 


.2  The  Maxwel 1 -Garnett 
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• 1 
prediction  of  the  classical  theory  was  discovered  by  Tanner  in  1972. 

The  experimental  FIR  absorption  coefficient  for  various  kinds  of 

metallic  particles  was  orders  of  magnitude  larger  than  the  prediction. 

Also,  he  observed  some  possible  evidence  of  the  quantum  size  effect 

in  small  particles  which  seems  to  arise  from  the  energy  level  correlation 

calculated  by  Gor’kov  and  Eliashberg.^  But,  to  the  present  time,  both 

the  classical  and  the  quantum  mechanical  calculation  have  failed  to 

give  the  correct  explanation  of  the  anomalous  FIR  absorption  by  small 

particles  inspite  of  the  agreement  of  the  frequency  dependence  with  the 

experimental  results  and  no  clear  answer  has  been  given  yet. 

Rev i ew 

Experiment 

The  first  FIR  measurements  on  a '^^100  A size,  gas-evaporated 
small  metallic  particles  were  done  by  Tanner  et  al.^^  Preparation  by 
gas-evaporation  was  expected  to  be  a promising  method  to  make 
well-characterized  spherical  small  particles  with  controllable  parameters 
such  as  the  average  particle  size,  possibly  the  particle  size 
distribution,  and  the  thickness  of  the  oxide  layer  on  the  surface  of  the 
metal  particles.  The  oxide  was  introduced  to  guarantee  the  electrical 
isolation  among  the  particles.  For  the  FIR  study  of  small  particles, 
they  used  the  samples  which  contained  free-standing  smoke  (clusters  of 
well-isolated  small  particles)  clamped  between  two  polyethelene  pieces. 
Because  the  density  of  the  smoke  was  very  low  in  comparison  with  that 
of  the  bulk  state,  the  volume  fraction  of  the  metal  in  the  sample  to 
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vacuum  was  less  than  0.1.  The  measured  FIR  absorption  coefficients 

of  small  particles  were  proportional  to  the  square  of  the  frequency, 

but  the  proportionality  constant  was  anomalously  larger  than  the 

theoretical  results.  Also,  no  temperature  dependence  of  the  FIR 

absorption  was  found  up  to  20°  K. 

1 5 

Granqvist  et  al.  repeated  the  FIR  measurement  on  aluminum 
small  particles,  confirmed  the  anomalous  FIR  absorption,  and  found  the 
absorption  to  increase  with  particle  size.  But  they  found  no  evidence 
of  the  Gor’kov  and  Eliashberg^  level  correlation  effect  and  they 
demonstrated  the  significance  of  the  particle  size  distribution  in  the 
observation  of  the  effect. 

A systematic  improvement  on  the  sample  preparation  was  made  by 
Russell  et  al.^^  They  were  able  to  control  the  volume  concentration  of 
the  small  metal  particles  by  using  an  insulating  medium,  KCl , as  a host. 
They  found  out  the  significance  of  the  eddy  current  loss  in  large  size 
particles  (1  ym)  and  the  linear  dependence  on  the  volume  fraction  in  the 
absorption  coefficient. 

At  the  same  period,  Carr  et  al.^^  measured  the  FIR  absorption  for 
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relatively  smaller  sizes  of  particles  than  those  of  Russell  et  al., 

but  large  enough  to  ignore  the  quantum  size  effect.  They  obtained 

similar  results  to  the  previous  measurements  except  that  the  absorption 

coefficient  of  aluminum  small  particles  showed  a quadratic  dependence  on 

volume  fraction.  They  also  measured  the  FIR  absorption  in  AI2O2. 

1 8 

which  was  a possible  explanation  for  the  anomalous  absorption,  and 
found  that  the  contribution  of  the  metal  oxide  was  not  responsible 
for  the  anomalous  absorption.  Except  for  the  Al  particles,  their 


FIR  data  fit 
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a 


(1.8) 


where  a is  the  absorption  coefficient,  is  a constant  in  cm,  f is 

the  volume  fraction  of  the  metal,  and  v is  the  frequency  in  cm"^. 

1 9 

Recently,  Devaty  and  Sievers  reported  a conclusive  result 

about  the  FIR  problem.  They  studied  a different  small  particle  system, 

100  A Ag  particles  embedded  in  a gelatin  matrix,  in  which  they  were  able 

to  control  the  dispersion  of  the  particles  and  found  out  that  a 

clustered  sample  had  a much  larger  FIR  absorption  coefficient  than  a 

well-dispersed  sample.  But  the  nonclustered  small  particle  sample 

still  showed  a factor  of  100  times  larger  absorption  than  theory.  Also, 

they  confirmed  that  the  absorption  coefficient  was  independent  of  both 

20 

temperature  and  magnetic  field  (up  to  60  kG).  However,  Lee  et  al. 
observed  a weak  temperature  dependence  of  the  FIR  absorption  for  fairly 
large  silver  particles  (a  > 500  A)  embedded  in  a Teflon  matrix.  Also, 
the  authors  measured  the  FIR  absorption  coefficient  which  is  close  to 
the  classical  predictions.  These  results  suggest  that  this  size 
particle  no  longer  retains  the  effect  of  the  broken  symmetry. 

The  early  work  on  the  optical  properties  of  small  particles  was 
started  by  investigating  the  color  change  in  a colored  NaCl  from  yellow 
to  blue  as  the  sample  was  annealed  in  1905.  Since  then,  many  studies 
have  been  made  on  small  particles  (colloids)  in  a colored  alkali  halide 
crystal.  The  interpretation  of  the  color  (F)  centers  and  the  colloids 
in  the  system  as  a one-component,  two-phase  thermodynamic  system  made  a 
successful  explanation  for  the  observed  phenomena.  In  other  words,  the 
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colloids,  precipitates  of  excess  alkali  metal,  behave  like  a condensed 

phase  and  the  F centers  behave  like  a vapor  phase.  The  blue  shift  in 

the  resonance  from  the  F center  peak  was  interpreted  as  the  corresponding 

Mie  resonance  (surface  plasmon  mode)  of  the  particles.  Another  system 

22 

used  for  the  study  of  the  surface  plasmon  was  Au-Si02  composite.  The 

gold  small  particles  were  grown  as  precipitates  by  irradiation  in  Si02 

matrix.  Although  in  general  there  was  a lack  of  information  about  the 

size  distribution  and  crystallinity  of  the  particles  grown  in  the 

matrix,  the  Mie  theory  was  qualitatively  successful,  even  though  the 

experiments  showed  more  broadening  in  the  resonance  curve  than  the 

, 22-31 

theory  predicted,  as  summarized  in  Table  1.1. 

Theory 

Classical  Theory 

We  want  to  calculate  the  electromagnetic  response  of  a particle 

and  the  cooperative  behavior  of  a collection  of  particles  dispersed  in  a 

medium.  The  response  of  a single  particle  can  be  calculated  by  the  Mie 

theory  or  by  the  approximation  of  the  quasi-static  limit.  The  details 

32 

of  the  Mie  theory  can  be  found  in  a review  paper  by  Carr  et  al . and 

the  references  therein.  However,  the  numerical  calculation  done  by 
33 

Henry  showed  that  the  difference  between  the  Mie  theory  and  the 

quasi-static  limit  is  negligible  for  the  small  size  particles  in  which 

we  are  interested  (a  100  A)  and  that  only  dipole  terms  are  important 

in  the  Mie  series.  The  quasi-static  limit  calculation  of  the  electric 

and  magnetic  polarizability  of  a small  particle  in  a dielectric  medium 

34 

can  be  found  either  in  Jackson  or  in  Landau  and  Lifshitz. 
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Table  1.1  Summary  of  optical  data 


Workers 


<a> 

(A) 


Av 


1 

r^Fi 

2ttc 

. 9 . 

^''obs 

^fit 

cm”^ 

) 

(cm“^ ) 

a) 

System 

Studied 


Scott  et 


al.23 

50 

912 

2,060 

22 

K-KBr 

Doyle^^ 

— 

161* 

1 ,610-4,840 

36-12 

Na-NaCl 

00 

Doremus 

26 

2,836 

2,760 

27 

Au-Glass 

60 

1,229 

2,140 

34 

Au-Glass 

100 

737 

3,560 

21 

Au-Glass 

Jain  and 

Arora^S 

75 

983 

3,810 

19 

Ag-KCl 

Kreibig^^ 

28 

2,634 

3,730 

20 

Au-Glass 

Granqvist 

and  Hunderi^/ 

16 

4,610 

8,440 

8.7 

Au  deposit 
on  substrate 

Norrman  et 

al.28 

^^50 

1,470 

5,680 

13 

Discontinuous 
Au  film 

Persson  and 

Liebsch29 

^-50 

1,475 

12,100 

6.1 

Ag-gelatin 

Gibson  and 

Buhrman^O 

^.50 

1,470 

14,300 

5.1 

Au-MgO 

granular  film 

'v50 

1,470 

7,580 

9.7 

Au-Al 2O3 
granular  film 

Kaiser  et 

al.31 

'v32 

2,304 

8,066 

9.1 

Discontinuous 
Au  film 

Note:  <a>:  average  particle  radius 

*•  ® ^bulk  " ^ 


Av: 

Av 


expected  broadening  at  FWHM 


i 


obs' 

fit- 


observed  broadening  at  FWHM 
electron  mean  free  path  used  to  fit  to  data 
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In  the  quasi-static  limit,  the  depolarization  field  can  be 
regarded  as  uniform,  determined  by  the  depolarization  factor.  For  a 
sphere,  we  can  write  the  electric  field  inside,  1^.^,  as 


(1.9) 


where  t is  the  applied  field  and  is  the  electric  polarizability  per 
unit  volume  of  the  sphere.  Also,  by  solving  the  standard  boundary 
value  problem,  we  get 


3e, 


'm 


+ 2e, 


t - 


m 


+ 2e^ 
m 0 


(1.10) 


where  e^,  are  the  dielectric  constant  of  the  medium  and  the  particle, 
respectively.  Therefore,  from  Equations  (1.9)  and  (1.10)  we  can 
identify  the  as 


Y 


e 


3_ 

4tt 


+ C- 


+ 2e, 


(1.11) 


One  important  remark  about  the  polarizability  y^  in  Equation 
(1.11)  is  that  y^  is  singular  when  e^(w)  = -2e^  which  predicts 
the  existence  of  the  surface  plasmon.  For  a Drude  metal  at  the  high 
frequencies. 


e^(co)  1 


2 


(1.12) 
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where  cOp  is  the  plasma  frequency  defined  as 
surface  plasmon  frequency  w^p  is 


f4„ne2lV2 


m 


Hence,  the 


Cl) 


sp 


(1.13) 


The  contribution  of  the  magnetic  dipole  moment  in  the  FIR  region 

is  caused  by  a non-zero  magnetization  in  the  medium  because  of  the  eddy 
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current.  This  contribution  considered  by  many  workers  » ’ ’ is 

known  to  be  important  for  a small  metal  sphere.  By  solving  the 
boundary  value  problem  of  a conducting  sphere  in  a time  varying  uniform 
magnetic  field,  the  magnetic  polarizability  per  unit  volume 
is 


Y 


m 


- -j^  cot  (ka) 


(1.14) 


where  k = /e~  - is  the  wave  vector  inside  of  the  particle. 

V m c 

In  order  to  describe  a system  with  small  metal  particles  in  a 

medium,  some  kind  of  averaging  process  is  necessary.  One  model  is  the 

Maxwell -Garnett  theory  (MGT),  an  analogy  to  the  Clausius-Mossoti 
37 

formula.  The  MGT  gives  the  effective  dielectric  function  of  the 
system  as 


^MGT  ~ ^0 
MGT  0 


'm 


e_ 

m 0 


(1.15) 
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where  is  the  effective  dielectric  function  and  f is  the  volume 
fraction  of  metal.  For  a sufficiently  dilute  system  (f  « 1),  Equation 
(1.15)  can  be  approximated  as 


'MGT  _ t ^ . ^16  2.2  2 , 

——  = 1 + 47Tf  Yg  + ^ TT  f Yg  + 0(f  ) 


(1.16) 


Similarly,  the  effective  magnetic  permeability  can  be  expanded  within 
MGT  as 


^MGT 


= 1 + 47Tf  Yj,  + O(f^) 


(1.17) 


Here  the  magnetic  permeability  of  the  medium  is  chosen  to  be  unity. 
Therefore,  the  absorption  coefficient  a|^^g-|-  is  obtained  from 


“MGT  ■ ^c^  /^MGtI^MGT 


(1.18) 
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Another  averaging  process  is  known  as  the  Bruggeman  theory  or 
the  Effective  Medium  Approximation  (EMA).^^’^^  The  EMA  asserts  that  the 
position-dependent  electric  field  in  the  medium  becomes 


t(r) 


c(r)  - e 


EMA 


c(r)  + 2c 


EMA 


< E> 


(1.19) 


therefore 


e - e 


< 


EMA 


e + 2c 


->  = 0 


EMA 


(1.20) 
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Here  e(r)  is  the  position-dependent  dielectric  function,  is  the 
effective  dielectric  function,  and  the  angular  bracket  represents  the 
ensemble  average.  Now,  if  we  assume  that  the  system  satisfies 
ergodicity,  the  condition  to  be  satisfied  to  have  a meaning  of  average 
for  a truly  random  system,  then  we  can  replace  the  ensemble  average 
with  the  volume  average  such  as 


The  EMA  thus  predicts  the  same  result  as  MGT  to  the  first  order  in  f. 

In  the  limit  of  low  frequency  and  small  volume  fraction,  the 
problem  is  now  reduced  to  the  calculation  of  the  susceptibilities  of  a 
small  particle  so  that  from  Equations  (1.16),  (1.17),  and  (1.18)  we  have 


<t>=ft  + (1  - f)t 

m ' 0 


(1.21) 


where  is  the  electric  field  in  particle  and  t is  the  field  in 

m ~ r\ 


'm  0 


medium.  From  Equations  (1.19)  and  (1.20)  immediately  we  get 


It  is  remarkable  that  the  EMA  predicts  the  metal-insulator 


transition  at  f = -j.  In  contrast,  the  MGT  has  no  transition  point. 
Again,  in  the  dilute  limit 


^EMA 


(1.23) 


c 


(1.24) 
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where  ^ » absorption  by  the  medium.  In  the  low-frequency 

limit,  the  imaginary  part  of  Vg  can  be  written  as 


ImYg 


9e  wa, 
0 I 

If;  2,  2 ^ 2x 

Idtt  (a-j  + 02) 


(1.25) 


with  e (w)  = 1 + and  a(w)  = a,  (co)  + io9(aj).  Also,  by  expanding 

the.equation  (1.14)  in  the  ka  « 1 limit  (however,  a little  care  is 
necessary  during  the  expansion^^ ) , we  get 


Imy. 


m 


40it 


Im(ka)‘ 


coa 


10c 


? ^1 


(1.26) 


Therefore,  after  subtracting  or  ignoring  the  absorption  by  the  medium, 
the  FIR  absorption  coefficient  in  MGT  becomes 

(1.27) 
or 


^ (jj  / — 

°^M6T  ^ 2/^0 

c 


9ca,e^ 
1 0 


4Tr(a^  + o^) 


. 2tt  ,2 
* 5?"  °1 


“MGT  = 


(1.28a) 


in  experimentally  more  appealing  form.  Here  is  the  proportionality 
constant  of  the  electric  dipole  term. 


K 


e 


(Oi/SOc) 

10[a-,/30c)^  + (02/300)^] 


(1.28b) 
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and  is  the  constant  for  the  magnetic  dipole  contribution, 

= 48TT^a^ey^  (a-,/30c)  (1.28c) 

with  V = (jo/27rc  and  a-j/SOc  in  fi'^cm"^. 

A test  of  whether  the  measured  FIR  absorption  follows  an 
3/2  1/2 

Cq(w)  ' or  an  e^(a))  ' dependence  can  be  done  to  identify  the 
importance  of  the  electric  dipole  and  the  magnetic  dipole  contributions 
to  the  absorption. 

Another  remark  is  that  there  is  an  enhancement  of  the  static 
dielectric  constant  of  the  small  metal  particle  composite.  From 
Equations  (1.16)  and  (1.23),  we  have 

Cgff  (1  + 3f)cQ  for  and  f « 1 (1.29) 


where  is  the  effective  dielectric  function  (either  or 
etr  ' EMA 

^MGT^’  ® different  point  of  view.  Equation  (1.29)  can  be 

understood  as  a result  of  the  expansion  of 


^eff 


(1.30) 


with  V = 1 and  f « f^  = 1/3.  Equation  (1.30)  shows  the  scaling  law 

behavior  near  the  critical  volume  fraction  f^.  The  actual  value  of  v 

and  f^  in  a three-dimensional  small-particle  composite  has  been 

measured  by  Grannan  et  al.^^  and  calculated  by  Straley^^  and  Bergman 
44 

and  Imry. 
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Quantum  Mechanical  Theory 

Since  the  kinetic  energy  of  electrons  is  quantized  as  discussed 
earlier,  the  standard  quantum  mechanical  calculation  of  the  induced 
dipole  moment  of  the  system  can  be  done  by  treating  the  interaction 

-V 

= eE  • r as  a perturbation.  Here  r is  the  position  vector  of  an 
electron  measured  with  an  origin  at  the  center  of  the  system,  ? is  the 
applied  electric  field,  and  e is  the  electron  charge.  Also,  the 
distinction  between  the  highly  degenerate  and  the  nondegenerate  energy 
level  contribution  to  the  dipole  moment  can  be  realized  as  follows.  To 
the  first  order,  the  perturbed  one-electron  wave  function  becomes 


+ y with!  = Ek  (1.31) 

where  is  the  unperturbed  state  eigenfunction,  E is  the 
corresponding  eigenvalue,  and  k is  the  z-direction  unit  vector. 
Therefore,  the  dipole  moment  p can  be  calculated  by 


Cfjezl'fV  = H 


2e' 


I 

m/n 


z n> 


^m- 


E. 


(1.32) 


In  order  to  visualize  the  effect  of  the  degeneracy,  let  us  consider  the 
excitation  between  adjacent  levels  with  average  energy  gap,  AE,  for 
simplicity.  The  the  dipole  moment  becomes 

p 'v  ^ |<1Iz|2>|2  E 


(1.33) 


24 


and  corresponding  static  electric  polarizability  is  obtained  by 
using  the  definition  p = Thus,  we  obtain 


a 

e 


=21 

aAE 


2 

a with 


l<l|z 


2>1^  a^ 


(1.34) 


Hence,  the  enhancement  on  the  electric  polarizability  due  to  the  level 
quantization  is 


aAE 


(1.35) 


where  is  the  classical  electric  polarizability,  a^  for  a metal 
sphere.  Now,  for  the  Kubo  system,  on  the  one  hand,  AE  is  the  order  of 
Ep/Ngi  therefore. 


a 


a 


2 N 


(1.36) 


with  N„ 
e 


’oJ 


, where  a^  is  the  Bohr  radius.  Consequently,  we  have  the 


enhancement  which  increases  quadratically  with  the  particle  radius  a. 
On  the  other  hand, for  a highly  degenerate  system  which  has  an  energy 
gap  order  of  Ep/N^^,  the  ratio  of  the  polarizability  gives 


0 mV3 
a 2 N„ 
e . e e 

_ % = 'Xj 

cl  a Er 

a r 

e 


e'/a. 


'V  1 


(1.37) 


which  shows  no  enhancement. 
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However,  in  a realistic  system,  the  one  electron  picture  breaks 
45 

down  because  of  screening.  In  the  classical  approach,  the  screening 
effect  can  be  included  by  replacing  the  external  electric  field  t by 
the  field  inside  of  the  particle  t.  as  in  the  following: 

^ = = (’-38) 


or 


%bs 


4tt 

3V  “e 


(1.39a) 


(1.39b) 


Here  is  the  observed  electric  polarizability,  Xg  is  the  electric 
susceptibility  defined  by  ^ = Xg^-jp*  ^nd  is  defined  as  the  dipole 

moment  per  unit  volume.  In  the  classical  limit  (AE  ->  0), 

cl  3 

approaches  the  classical  value  (a^  = a ) which  satisfies  the 
correspondence  condition. 

An  alternative  approach  to  the  screening  can  be  made  by 

considering  the  Coulomb  interactions  between  electrons  at  different 

energy  levels,  even  though  the  screening  in  a small  particle  is  smaller 

than  that  in  the  bulk  state.  Lushnikov  and  Simonov^^  included  the 

electron-electron  interaction  in  the  calculation  of  a within  the 

e 

Random  Phase  Approximation  and  found  out  that  1 for  a small 

particle  at  metal  density.  These  results  explain  the  reason  for  the 
absence  of  the  enhancement  on  the  observed  electric  polarizability. 
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But,  there  still  exists  a possibility  of  experimental 
observation  of  the  effect  of  the  energy  level  correlation  which  was 
first  realized  by  Gor’kov  and  Eliashberg.^  In  a recent  calculation, 
Devaty  and  Sievers^^  showed  that  a nonquadratic  frequency  dependence  on 
the  FIR  absorption  coefficient  could  be  regarded  as  the  Gor’kov  and 
Eliashberg  effect.  However,  the  magnitude  of  theoretical  FIR 
absorption  coefficient  is  still  orders  of  magnitude  smaller  than  the 
experimental  results. 

40 

Kawabata  and  Kubo  used  the  fluctuation-dissipation  theorem  to 
calculate  the  optical  conductivity  of  a small  particle  in  visible 
frequency  range  and  gave  a quantum  mechanical  interpretation  of  the 
broadening  of  the  surface  plasmon  resonance  peak.  Their  calculation 
showed  that  the  optical  conductivity  of  a small  particle  was  broadened 

significantly  from  the  Drude  conductivity. 

4Q 

Wood  and  Ashcroft  also  carried  out  a calculation  of  the 
dielectric  function  of  a small  metallic  cube  by  evaluating  the 
momentum  matrix  element  explicitly  within  the  linear  response  theory 
and  the  random  phase  approximation.  They  obtained  a similar  result  as 
Kawabata  and  Kubo  except  for  a geometrical  factor.  However,  their 
result  failed  to  explain  the  anomalous  FIR  absorption. 

The  small  particle  problem  can  be  summarized  as  follows: 

1.  The  FIR  absorption  coefficient  can  be  expressed  as 
2 

a = Kg^p  fv  except  for  the  aluminum  small  particle  case. 
*^exp  *^theory. 
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3.  K is  not  a function  of  either  temperature  or  magnetic 

exp 

field. 

4.  K is  increased  as  the  particle  size  increases. 

Search  for  the  Solution 

The  following  questions  are  raised  in  order  to  understand  the 
anomalous  FIR  response  of  the  small  particle  system  from  experimental 
and  theoretical  points  of  view:  whether  the  anomalous  FIR  absorption 

wil-1  be  persistent  to  mid-infrared  (MIR)  frequencies  or  not,  whether 
the  clustering  is  responsible  for  the  observed  anomalous  FIR  absorption 
or  not,  whether  there  exists  a possible  absorption  mechanism 
incorporated  with  phonons  or  not,  and  whether  the  choice  of  the  Drude 
type  dielectric  function  is  adequate  for  a small  particle  or  not. 

Clustering  Effects 

Many  efforts  have  been  made  to  understand  the  anomalous  FIR 

absorption  by  small  particles.  One  possible  candidate  is  clustering. 

Because  there  was  no  convincing  direct  experimental  evidence  on  the 

randomness  of  the  dispersed  small  particles  inamedium,  the  effect  of 

small  particle  clusters  could  play  an  important  role  in  the  FIR 

absorption  mechanism. 

50-53 

Many  workers  have  demonstarated  that  the  FIR  absorption  can 

be  enhanced  by  clustering  to  be  close  to  the  experimental  results. 

This  enhancement  occurs  because  of  the  poor  effective  conductivity  of 

individual  clusters.  However,  in  the  clustering  model,  the  quadratic 

frequency  dependence  inevitably  breaks  down  as  shown  by  Sen  and 
54 


Tanner. 
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In  this  work,  more  systematic  analysis  on  sample  characterization 
and  FIR  measurements  will  be  made  for  clustered  samples  and  presumably 
nonclustered  samples. 

Phonons 


in  a small  particle  by  direct  coupling  of  external  electric  field  to 
the  phonon  through  the  unscreened  surface  ions.  The  theory  was  able  to 
show  the  temperature  independence  and  enhanced  absorption  but  gave  a 
wrong  particle  size  dependence  on  the  absorption  coefficient.  Also, 
an  effect  of  surface,  which  modifies  the  phonon  density  of  states,  was 
not  included. 

Another  possibility  is  the  Holstein  process  in  a small  particle, 

since  the  particle  dimension  is  much  less  than  the  classical  skin 

depth  and  electron-phonon  coupling  gives  a frequency-dependent 

relaxation  time  in  addition  to  the  lifetime  in  the  bulk  state.  For  a 

small  particle,  the  phonon  density  of  states  has  to  be  modified  by 

including  the  surface  contribution.  By  solving  the  Weyl's  boundary 
56 

value  problem,  we  obtain  the  phonon  density  of  states  such  as 


55 

Click  and  Yorke  suggested  a possible  enhanced  FIR  absorption 


(1.40) 


where 
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Here  and  are  the  transverse  and  the  longitudinal  sound  velocity, 

respectively.  With  the  modified  density  of  states  and  renormalized 

plasma  frequency,  we  can  evaluate  the  electron-phonon  contribution  in 

57 

FIR  by  making  use  of  Allen's  formula.  However,  the  results  show  a 

58 

quartic  dependence  on  frequency  in  the  absorption  coefficient,  which 
is  not  consistent  with  the  observation.  It  is  likely  that  the 
electron-phonon  interaction  would  be  less  effective  in  a small  particle 
because  of  the  difficulties  in  energy  and  momentum  conservation  due  to 
the  quantization  of  electron  energy  and  phonon  frequency. 

Non-Drudian  Possibility 

The  last,  conceptually  rather  difficult,  possibility  is  the 

59  • 

non-Drudian  behavior.  For  example,  the  Anderson  localization  is 
generally  believed  to  occur  in  an  electron  system  with  impurities  whose 
concentration  is  above  a certain  critical  concentration.  We  can 
interpret  the  diffusionless  process  as  a consequence  of  locally  broken 
symmetry  in  the  translational  motion  of  electrons. 

Now,  for  a small  particle,  the  quantum  mechanical  role  of  the 
surface  localizes  the  electrons  and  quantizes  their  kinetic  energy.  In 
fact,  the  finiteness  of  the  size  (or  broken  symmetry)  imposes  another 
important  consequence;  unlike  the  case  of  the  electron  gas  in  the  bulk 
state,  the  momentum  matrix  element  does  not  vanish.  To  put  this  idea 
another  way,  the  position  matrix  element  of  an  electron  is  always  finite 
in  such  a finite  system,  requiring  a fundamental  modification  to  the 
calculation  of  the  conductivity  of  a small  particle. 


30 


These  points  will  be  extensively  studied  experimentally  and 
theoretically  in  this  dissertation. 

In  the  next  chapter,  detailed  experimental  techniques  will  be 
discussed  including  sample  preparation  and  characterization, 
instrumentation,  and  measurements. 

In  Chapter  III,  the  experimental  measurements  of  FIR,  MIR,  and 
visible  frequency  response  of  Al-KCl,  Au-KCl,  and  Ag-KCl  small  particle 
composites  and  the  interpretation  of  the  results  will  be  given. 

A quantum  mechanical  theory  of  electrical  and  optical 
conductivity  of  a small  particle  will  be  presented  along  with  the 
energy  level  correlation  effect  and  a possible  explanation  for  the 
broadening  of  the  surface  plasmon  mode  will  be  made  in  Chapter  IV. 

Finally,  a summary  and  conclusion  will  be  given  in  Chapter  V. 


CHAPTER  II 
EXPERIMENT 


Sample  Preparation  and  Characterization 

Small  Particles 

Small  particles  were  prepared  by  the  gas-evaporation  technique. 

As  shown  in  Figure  2.1,  small  particles  were  deposited  on  the  inside  of 
a fused  quartz  cylinder  (5-1/2"  I.D.  and  12"  long)  in  the  form  of  smoke 
(aggregates  of  small  particles)  by  the  evaporation  of  a metal  (A1 , Ag, 
or  Au)  under  an  inert  gas  atmosphere. 

There  are  three  major  controllable  factors  which  determine  the 
size  of  the  particle  and  the  size  distribution  in  this  method: 

1.  Type  of  inert  gas 

2.  Pressure  of  the  inert  gas  in  the  vacuum  chamber 

3.  Evaporation  rate  (temperature  of  the  source) 

The  first  and  the  second  condition  mainly  determine  the  size  of  particles. 
The  average  size  of  particles  is  proportional  to  the  pressure  and  the 
atomic  mass  of  the  inert  gas  used  and  the  vapor  pressure  and  the  atomic 
mass  of  the  metal.  Throughout  this  experiment.  He  or  Ar  gas  was  used  for 
A1  small  particles  and  Ar  + 10%  O2  gas  mixture  was  used  for  Ag  and  Au 
small  particles.  The  purpose  of  the  oxygen  was  to  guarantee  the 
electrical  isolation  between  particles  by  forming  an  oxide  layer  on  the 
surface  of  the  particles.  The  background  pressure  before  introducing 
the  inert  gas  was  less  than  5 x 10”^  Torr.  The  gas  pressure  of  the 
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Figure  2.1  Small  particle  preparation  apparatus 
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inside  of  the  bell  jar  was  controlled  by  changing  the  flow  rate  and 
the  pumping  speed.  The  third  condition  mainly  contributes  to  the 
yields  and  the  size  distribution.  As  a heat  source,  several  different 
designs  of  tungsten  baskets  were  sued.  Among  them,  a basket  with  two 
tungsten  wires  (9.5  mm  diameter)  twisted  together  gave  the  best  results. 
Also,  an  Alumina-coated  Mo  boat  was  tested  for  comparison.  Typically, 
about  100  200  mg  of  pure  A1  (99.9999%),  Ag  (99.9995%),  or  Au  (99.99%) 

was  used  to  prepare  small  particles. 

A Transmission  Electron  Microscope  (TEM)  was  used  to  measure  the 
particle  size  and  the  size  distribution.  Also,  the  thickness  of  the 
oxide  layer  could  be  measured  in  principle  by  comparing  the  bright  field 
picture  of  a particle  with  the  corresponding  dark  field  picture  to  ± 3 A 
accuracy.  The  typical  thickness  which  is  known  to  be  the  order  of  10  A 
within  the  instrumental  resolution  (3  A line  to  line).  In  addition,  the 
electron  diffraction  and  the  EDX  (Energy  Dispersive  X-ray)  spectroscopy 
analysis  of  small  particles  were  done  through  the  TEM.  The  TEM  pictures 
of  two  different  sizes  of  A1  small  particles  are  shown  in  Figures  2.2 
and  2.3  and  a result  of  the  small  particle  data  is  given  in  Table  2.1. 
Also,  the  result  of  the  oxide  thickness  estimation  is  shown  in  Table  2.2. 

Composites 

Pure  KCl  was  used  as  a dielectric  medium  for  all  frequency 
ranges  (FIR,  MIR,  Near-IR,  and  Visible)  measurements.  Also,  paraffin 
was  used  for  FIR  measurements  over  the  frequencies  of  the  optical  phonon 
band  of  KCl . 
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Figure  2.2  Transmission  electron  micrograph  of  A1  15  (top)  and  the 
corresponding  electron  diffraction  pattern  (bottom) 


Arb.  Scale 
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0-10  KEV 

Fiijure  2.3  Transmission  electron  micrograph  of  A1  18  (top)  and  the 
corresponding  EDX  analysis  (bottom)  which  shows  no 
contdiiiinaticn  by  foreign  atuiiiS  within  the  resolution 
(1  atomic  %) 


Table  2.1  Small  particle  parameters 
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Table  2. 1— Continued 
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Table  2.2  Analysis  of  the  oxide  layer 


Aluminum 

Particle  radius 
in  bright  field 

re(«) 

Particle  radius 
in  dark^field 
r„(A) 

Oxide  thickness 
t(X)  i iH 

(t  = rg-rp) 

A1  1 

100 

85 

15 

135 

120 

15 

80 

65 

15 

A1  9 

85 

75 

10 

95 

85 

10 

50 

40 

10 

50 

40 

10 

A1  13 

115 

100 

15 

70 

55 

15 

130 

115 

15 
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The  bulk  KCl  or  paraffin  was  ground  up  to  1 y size  grains  in  a 
Freezer  Mill  operating  at  77°K  and  then  the  dried  KCl  powder  (dried  at 
150°C  for  3 5 hours  under  vacuum)  or  paraffin  (dried  under  vacuum  at 
room  temperature)  was  mixed  with  small  particles.  The  volume  fraction 
(f)  was  estimated  by  (to  the  first  order) 


where  (p^)  is  the  density  of  metal  (host  medium)  and  (W^)  is  the 
weight  of  metal  (host  medium).  The  mixture  was  mixed  in  a Freezer  Mill 
and  then  pressed  in  an  evacuated  die  into  a 5/8"  diameter  pellet.  The 
compression  took  place  at  about  130,000  psi  for  KCl  and  at  3,200  psi  for 
paraffin.  The  composite  pellet  was  reground  several  times  until  a 
satisfactory  homogeneity  was  achieved.  Normally  four  times  of  regrinding 
was  enough  for  low  volume  fraction  samples  (f  % 0.008).  For  the 
visible  frequency  measurements,  the  samples  were  diluted  ten  times  from 
f = 0.001  samples  to  have  enough  transmission  of  light.  The 
pressure-regrinding  process  was  repeated  four  times  for  these  diluted 
samples.  The  thickness  of  the  sample  ranged  from  0.18  mm  to  5.0  mm, 
depending  on  the  volume  fraction  and  the  frequency  ranges  covered. 

Samples,  for  two  different  A1  particle  sizes  which  were 
believed  to  have  large  clusters  of  small  particles,  were  deliberately 
prepared  without  using  the  Freezer  Mill  in  order  to  study  the  effect  of 
clustering. 
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A correction  of  the  volume  fraction  is  necessary  to  estimate  the 
true  metal  volume  fraction  because  each  metal  particle  has  about  10  K 
thick  oxide  coating.  In  actual  calculation,  the  following  correction 
can  be  made  to  the  apparent  volume  fraction 


where  f is  the  corrected  volume  fraction,  a is  the  particle  radius,  and 
t is  the  thickness  of  the  oxide  coating. 

Also,  in  order  to  get  rid  of  the  interference  fringes  in  the 
FIR  transmission  measurement  due  to  the  multiple  internal  reflection 
between  front  and  back  surfaces  of  a sample,  one  series  of  samples  had 
an  extra  layer  of  KCl  added  to  make  the  total  thickness  of  the  sample 
larger  than  about  5 mm.  Under  the  assumption  that  the  difference  between 
the  static  dielectric  constants  of  the  small  particle-KCl  composite 
and  of  the  pure  KCl  layer  is  negligible,  the  FIR  beam  will  not 
distinguish  the  boundary. 

Sample  Characterization 

A TEM  or  an  optical  microscope  was  used  to  characterize  the 
randomness  of  small  particle  distribution  in  a composite.  For  the  TEM 
study,  very  thin  sections  of  a sample  were  made  by  using  either  an  ion 
mill  or  a chemical  polisher.  In  order  to  mill  down  a sample  through 
the  ion  mill,  a sample  was  cut  into  3 mm  diameter  pellets  and  ground 
down  to  100  y thickness  and  dimpled  until  the  thickness  of  the  thinnest 
portion  became  about  10  y.  Then,  the  ion  mill  was  used  to  make  a pinhole 
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on  the  specimen.  The  TEM  pictures  were  taken  in  the  vicinity  of  the 
pinhole  where  the  electrons  can  be  transmitted. 

Another  effective  method  was  attempted  by  using  a chemical 
polisher.  Since  water  dissolves  KCl , a 90%  alcohol  and  10%  distilled 
water  mixture  was  used  as  a chemical  solution.  A 3 mm  diameter  sample 
was  sandwiched  between  two  jet  streams  of  alcohol -water  mixture  and 
polished  down  until  a pinhole  was  made. 

It  took  about  3-5  minutes  to  prepare  a sample  with  the  chemical 
polishing  method  in  contrast  to  5-10  hours  in  the  ion  milling  method. 

The  prepared  sample  was  placed  on  a cold  stage  operating 
normally  at  -100°C  to  prevent  the  melting  of  the  sample  when  the 
electron  beam  was  condensed  on  a particular  area,  and  the  sample  stage 
was  mounted  in  the  TEM.  About  10-15  different  areas  were  examined  for 
each  sample  and  the  results  are  shown  in  Figures  2.4,  2.5,  and  2.6. 
Individual  particles  can  be  seen  clearly  in  the  TEM  pictures. 

Especially  because  the  typical  thickness  of  the  area  examined  is  about 
1,000  A,  the  composite  is  surprisingly  close  to  the  truly  random  system 
with  well  separated  individual  particles  for  samples  made  with  A1  14 
(103  A)  particles. 

In  addition,  an  optical  microscope  and  TEM  study  was  done  on  the 
clustered  samples  which  were  mixed  without  using  the  Freezer  Mill  to  see 
the  size  and  the  distribution  of  clusters  as  shown  in  Figures  2.7  and 

2.8. 
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Figure  2.4  Transmission  electron  micrographs  AI-KCl  composites 

(a)  Water  polished  f = 0.032  A1  1-  Kcl  composite  plate  (# 
of  regrinding  = b)  and  (b)  ion  milled  f = o.OOl  A1  18-KCl 
composite  plate  with  EDX  result  which  shows  some  Fe 
contunination  (circled  area)  during  the  sample  manufacturinn 
process  {%  of  regrinding  = 4) 
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2000A 

Figure  2.5  Optical  and  trunsniission  electron  micrographs  of  f = 0.001 
A1  18  - KCl  composite  (#  of  regrinding  = 4) 

(a)  Optical  micrograph  of  chemically  polished  plato  and 

(b)  transmission  electron  micrograph  of  ion  milled  plate 
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Figure  2.6  Transmission  electron  micrugrapiis  uf  A1  14-KCl  composite 
(f  = 0.001,  ground  4 times) 

(a)  Ion  milled  and  (b)  chemically  polished  plate 
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(b)  r 1 
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Figure  .?.7  Optical  and  transmiss ion  electron  micrographs  of  A1  18-KCl 
composite  (f  = 0.001,  mixed  by  hand) 

(a)  Optical  microscope  result  and  (b)  TEM  result 
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rigure  2.8  Optical  anu  ti unsmissicn  electron  micrographs  of  A1  14-KCl 
composite  (mixed  by  hand) 

(a)  .Optical  microscope  result  (f  = 0.002)  ^nd  (b)  TEM 
result  (f  = 0.001 ) 
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Infrared  and  Optical  Technique 
Fourier  Transform  Spectroscopy 

Interference  of  light  from  two  mutually  coherent  sources  is  a 
well-known  phenomenon.  As  illustrated  in  Figure  2.9(a),  the  intensity 
I at  D can  be  written  as  the  time  average  of  the  square  modulus  of  the 
sum  of  two  amplitudes  A-|(t)  and  A2(t  + t),  which  correspond  to  the 
amplitude  from  slit  1 and  slit  2,  respectively. 

I = <[A^(t)  + A2(t  + T)][A^(t)  + A2(t  + t)]*> 

= I^  + l2  + 2Re  <A^(t)  A^(t  + x)>  (2.3) 

where  I^  = <lA-|(t)l^,  I2  = <|A2(t  + t)|^,  and  the  angular  bracket 
represents  the  time  average.  I^  (I2)  is  the  intensity  of  light  through 
SI  (S2)  only.  The  last  term  in  Equation  (2.3)  contains  all 
interference  effects  because  of  the  temporal  correlation  between  two 
amplitudes  caused  by  t.  When  t “ (i.e.,  h ->  “),  there  is  no 
correlation.  Therefore, 

1=1  = I,  + I5  (2.4) 

Now,  if  we  can  manage  to  give  the  time  difference  t between  two  light 
beams  in  the  limit  of  h 0,  then  |A-|  1 = IA2I  by  ignoring  the 
fluctuation  in  the  amplitude.  Hence, 


I = I„  + 2Re  <A(t)  A*(t  + t)> 


(2.5) 
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(a) 


Figure  2.9  Two  different  illustrations  of  obtaining  interference 
(a)  Double  slit  and  (b)  Michel  son  interferometer 
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where  the  subscripts  1 and  2 have  been  deleted.  The  function 

<A(t)  A*(t  + t)>  is  called  as  the  autocorrelation  function;  it  satisfies 

the  Wiener-Khintchine  theorem. That  is, 

<A(t)  A*(t  + t)>  = FT  [ A(v)  (2.6) 

where  FT  represents  the  Fourier  transformation  and  v represents 
frequency.  As  shown  in  Figure  2.9(b),  the  Michelson  interferometer  is  a 
perfect  example  of  what  can  be  described  by  Equation  (2.5).  By 
changing  from  the  time  domain  to  a spatial  coordinate,  such  as  6 = ct 
where  6 is  the  path  difference  between  the  two  mirrors,  and  defining 
the  correlation  function  r(6)  as 

r(6)  = 2 Re  <A(0)  A*(6)>  (2.7) 

we  can  rewrite  expression  (2.5)  as 

1(6)  = + r(6)  (2.8) 

Again,  r(6)  satisfies  the  Wiener-Khintchine  theorem.  Therefore,  the 
power  spectrum  can  be  obtained  by  an  inverse  Fourier  transform  of  the 

r(6): 

00 

S(v)  = I r(6)  e^^^^'^d6  (2.9) 

-00 

r\ 

where  v = 2^  ^nd  S(v)  = |A(v)|  . Now,  when  6 = 0,  = 21^;  therefore, 

we  can  rewrite  Equation  (2.9)  in  the  desired  form 
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S(v)  = 


(1(6)  - 1 l^) 

-CO 


(2.10) 


1(6)  is  the  measured  quantity  in  the  experiment. 

In  actual  measurements  apodization  is  necessary  to  remove  the 

difficulty  arising  from  finite  maximum  path  difference.  Proper 

filtering  is  also  essential  to  overcome  the  aliasing  problem  due  to 

discrete  sampling  by  removing  the  frequencies  higher  than  the  maximum 
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cut-off  frequency  determined  by  the  sampling  interval. 


Optical  Spectroscopy 

As  the  frequency  range  of  interest  increases  towards  the  higher 
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frequency  region,  the  Fellget  advantage  which  counts  for  the  thermal 
detector  noise  only  in  the  infrared  loses  its  importance  due  to  the 
increasing  photon  noise  in  the  radiation  field.  Therefore,  a 
monochromator  is  normally  used  to  avoid  many  difficulties  in  Fourier 
transform  spectroscopy  in  the  Near-IR  and  visible  frequency  range. 

But,  the  photon  noise  also  eventually  limits  the  performance  of  the 
monochromator  at  very  high  frequencies. 

Generally,  a grating  monochromator  is  used  by  applying  the  rule 
of  diffraction.  For  a wave  length  X, 


= JL  1 
2d  sinG 


(2.11) 


where  n is  the  integer  order,  d is  the  grating  constant  (cm/line),  0 is 
the  angle  of  the  grating  measured  from  the  position  for  normal  incidence 
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of  the  light.  The  resolution  is  determined  by  the  slit  width  and  A0, 
which  is  the  angle  of  rotation  at  each  step. 

Instrumentation 

Michel  son  Interferometer 

In  the  FIR  region,  the  Fellget  advantage  is  obtained  from  the 
signal-to-noise  ratio  argument  and  Fourier  transform  spectroscopy 
maximizes  the  efficiency  of  the  experiment. 

To  measure  the  spectrum  in  the  10  to  800  cm”^  region,  a 
Michel  son  interferometer^^  was  used.  Correct  combination  of  the 
thickness  of  a Mylar  beam  splitter,  filters,  and  the  sampling  step  size 
is  important  to  remove  the  aliasing  problem.  As  shown  in  Figure  2.10, 
the  entire  instrument  is  placed  inside  a vacuum  tank  to  prevent  FIR 
absorption  due  to  water  vapor  in  the  air. 

Lamellar  Grating  Interferometer 
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A lamellar  grating  interferometer  was  used  to  cover  the 
frequency  range  from  2 cm'^  to  60  cm"^ . In  principle,  this  instrument 
works  in  the  same  way  as  the  Michel  son  type  interferometer,  but  has 
better  modulation  than  that  of  the  Michelson  interferometer,  which  has 
poor  beam  splitter  efficiency  in  the  very  FIR.  A detailed  instrumental 
set-up  is  shown  in  Figure  2.11. 

Grating  Monochromator 

A Perkin-Elmer  model  grating  monochromator  was  used  (Figure 
2.12)  to  cover  from  MIR  to  ultra-violet  frequency.  As  described  in  the 
optical  spectroscopy  section,  each  single  frequency  is  selected  by  using 
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Figure  2.10  Schematic  diagram  of  Michelson  interferometer 
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(b)  Cross  sectional  view  of  Lamellar  Grating 


Figure  2.11  Lamellar  grating  interferometer  (a)  and  cross  sectional 
of  the  lamellar  grating  to  show  the  optical  path 
difference  (b) 
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Figure  2.12  Schematic  diagram  of  grating  monochromator 
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the  grating,  filter,  and  the  matching  source.  A more  detailed  source 

of  information  about  the  operation  of  this  instrument  is  found  in 
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the  work  by  Cummings. 

Measurement 

Source  and  Detector 

A mercury  arc  lamp  was  used  as  a FIR  source,  and  the  source  was 
chopped  to  remove  all  background  radiation.  A doped  germanium  bolometer 
operating  at  1.2°K  served  as  a FIR  detector.  Figure  2.13  shows  a 
diagram  of  the  cryostat. 

In  order  to  detect  the  FIR  signal  from  the  instrument,  a 1/2" 

I.D.  brass  tube  was  used  to  guide  the  radiation  through  the  cryostat. 

For  higher  frequency  measurements,  three  different  types  of 
detectors  for  three  different  matching  sources  were  used  with  the 
Perkin-Elmer  grating  monochromator.  Table  2.3  shows  a summary  of  the 
combination  of  instrument,  source,  detector,  and  frequency  range  used 
in  this  work. 

Data  Acquisition 

The  FIR  signal  was  detected  by  a bolometer,  which  was  biased  from 
the  external  DC  voltage  source.  The  voltage  change  across  the  detector 
due  to  the  change  of  the  resistance  in  response  to  the  incoming 
radiation  was  fed  into  a look-in  amplifier  (Ithaco  393).  The  amplified 
signal  was  sent  to  a chart  recorder  and  a digital  voltmeter  (DVM).  The 
signal  was  averaged  by  the  DVM  (Fluke)  for  a given  time  interval 
(typically,  1 to  5 seconds,  which  depends  on  the  signal-to-noise 
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^ first  floor  sample  rotator 


Figure  2.13  Schematic  diagram  of  detector  cryostat 
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Table  2.3  Summary 


Frequency  (cm 


Source 


Detector  Instrument 


2-60 

Mercury  arc 
lamp 

Bolometer 

Lamellar  gratinc 
interferometer 

10-800 

Mercury  arc 
lamp 

Bolometer 

Michelson 

interferometer 

690-4,400 

Globar 

Thermocouple 

Perkin-Elmer 

monochromator 

4,400-15,400 

Tungsten 

lamp 

PbS 

Perkin-Elmer 

monochromator 

14,600-27,000 

Tungsten 

lamp 

Photomultipl ier 

Perkin-Elmer 

monochromator 

26,000-47,000 

D2  lamp 

Photomultiplier 

Perkin-Elmer 
monochroma tor 
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condition).  And  then,  the  digitized  information  by  the  DVM  was  stored 
in  a POP  11/23  computer.  After  having  taken  a point,  the  computer  sent 
a signal  to  the  stepping  motor  controller  to  advance  to  the  next 
step.  This  process  was  repeated  until  the  desired  number  of  data  points 
had  been  taken.  The  stored  data  were  analyzed  through  the  Fast  Fourier 
Transform  programs  after  having  finished  the  data  acquisition. 

For  the  grating  monochromator,  the  spectrum  was  directly 
acquired  from  the  detector  used  and  the  same  procedure  as  in  the  FIR 
measurement  case  was  repeated  until  a whole  spectrum  range  was  covered. 
The  spectrum  was  normalized  and  analyzed  through  the  computer. 


CHAPTER  III 
RESULTS 


The  electromagnetic  response  of  the  small  particle  composites 
were  analyzed  from  transmission  measurements,  with  the  absorption 
coefficient  given  by 


here  a(w)  is  the  absorption  coefficient  at  frequency  w,  d is  the 
thickness  of  the  sample,  T(w)  is  the  measured  transmission,  and  R(w)  is 
the  reflectance  of  the  sample.  Typically,  the  reflectance  of  a 
composite  sample  (f  = 0.03)  is  close  to  that  of  pure  KCl  and 
essentially  frequency  independent.  Therefore,  the  reflectance  of  pure 
KCl  was  used  in  the  data  analysis. 


Absorption  Coefficient 

Figure  3.1  shows  a result  of  the  FIR  measurements  on  small 
aluminum  particle-KCl  composites.  The  data  can  be  shown  to  fit  the 
following  formula: 


a(u)  = - T((jj)  + iln  [1  - R((jj)] 


(3.1) 


FIR  Data 


(3.2) 
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Figure  3.1  FIR  absorption  coefficients  of  A1  (A1  15)-KC1  composites 
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where  a is  the  correction  for  the  reflection  and  the  detector 
0 

nonlinearity.  The  proportionality  constant  can  be  a function  of 
a number  of  variables;  for  example, 


a,  f) 


(3.3) 


where  e (e  ) is  the  dielectric  constant  of  the  metal  (surrounding 
m 0 

medium),  a is  the  radius  of  the  particle,  and  f is  the  volume  fraction. 

As  shown  in  Figure  3.2,  the  dependence  is  quite  obvious  inthe 

a versus  plot.  In  this  particular  measurement,  an  extra  Kcl  layer 
exp 

was  used  to  remove  interference  effects,  as  described  earlier.  One 
remark  on  these  data  is  that  the  volume  fraction  dependence  of  the 
absorption  coefficient  is  close  to  linear  in  f in  contrast  to  the 
previous  work  done  by  Carr  et  al.^^  In  other  words,  K^^p  is  not  a 
function  of  f.  The  f^  dependence  in  the  previous  results  could  be 
attributed  to  the  effect  of  interparticle  tunneling  due  to  the 
clustering  of  particles. 

Another  experimental  characteristic  of  the  FIR  absorption  by 
small  particles  is  that  there  is  no  clear  distinction  among  different 
types  of  metals  which  suggests  that  the  FIR  absorption  mechanism  is 
governed  by  electrons.  For  example,  some  results  of  the  FIR 
measurements  on  Ag  and  Au  small  particles  show  about  the  same  order  in 
magnitude  of  absorption  as  that  of  the  similar  size  aluminum  particles 
and  also  show  the  quadratic  frequency  dependence  f the  absorption 
coefficient  as  shown  in  Figure  3.3. 


O' ( cm" ' ) 
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Figure  3.2 


ot0xp  vs.  V 
(A1  18)-KC1 


plot  of  FIR  absorption  coefficient  of  AL 
compos intes 
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Figure  3.3  FIR  absorption  coefficient  of  Ag-CKI  and  Au-KCl  composites 
(note  that  the  interference  fringes  are  arising  from  the 
multiple  reflection) 
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A comparison  of  the  experimental  absorption  coefficient  with 
the  result  of  the  corresponding  classical  theory  is  made  in  Figure  3.4. 
The  persistent  anomalous  FIR  absorption,  which  does  not  agree  with  the 
absorption  described  by  the  classical  small  particle  model  can  be 
clearly  seen  in  this  plot. 

The  particle  size  dependence  on  the  FIR  absorption  strength, 

K , in  Figure  3.5  shows  a difference  to  the  theory  which  suggests  that 
exp  ^ 

the  intrinsic  electrical  properties  of  a small  particle  could  be 
different  from  what  we  have  believed  to  be  correct  (e.g.,  the  concept 
of  small  "Drude  metal"  particle  and  the  surface  scattering  rate 
approximation,  etc.).  Also,  as  shown  in  Figure  3.6,  a correction  due 
to  the  oxide  coating  on  small  particles  was  made  on  under  the 
assumption  that  the  core  metal  part  of  particles  is  solely  responsible 
for  the  anomalous  FIR  absorption  by  calculating  the  true  volume 
fraction  such  as 


(3.4) 


where  a is  the  particle  radius  and  t is  the  thickness  of  the  oxide 
coating,  here  assumed  to  be  10  K. 


Clustering  Effect-FIR 

One  conjecture  made  to  resolve  the  FIR  problem  attributes  it  to 
clustering  of  the  individual  particles,  because  there  was  no  clear 
experimental  verification  of  randomness.  This  clustering  theory  does 
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ppQquGncy  (cm”*  ) 


Figure  3.4  A comparison  of  an  experimental  result  to  the  result  of 
classical  calculation 


K (cm) 
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Figure  3.5  Kexp  in  cm  vs.  particle  size  plot  along  with  the 

theoretical  calculation  of  Kg  (electric  dipole)  and 
(magnetic  dipole) 


(cm) 
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Figure  3.6  Kgxp  (closed  circle)  vs.  particle  size  plot;  also,  Kexp 
(opfen  circle)  is  shown 
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show  as  much  enhanced  FIR  absorption,  which  mainly  results  from  the 
poor  effective  conductivity  of  clusters  (order  of  micron  size)  as 
compared  to  a metallic  particle. 

From  the  experimental  point  of  view,  two  tests  of  the  mixing 
and  grinding  effect  on  the  small  particle-insulator  composites  can  be 
done:  first,  directly  by  TEM  analysis,  and,  second,  by  the  FIR  and  MIR 

measurements  on  these  samples.  As  shown  in  Chapter  II,  the  results  of 
the  TEM  analysis  showed  no  clustering  in  carefully  prepared  composites 
and  showed  a clear  difference  between  the  random  and  the  clustered 
samples.  The  effect  on  the  FIR  absorption  of  successive  regrinding 
by  the  "freezer  mill"  is  shown  in  Figure  3.7.  The  absorption 
coefficient  of  sample  mixed  by  hand  also  is  shown  in  Figure  3.7.  The 
FIR  absorption  strength,  K » decreases  significantly  after  first 

viA  [J 

regrinding  by  the  "freezer  mill."  However,  the  grinding  effect  shown 
in  Figure  3.7  demonstrates  a tendency  of  increasing  which 

saturates  beyond  a certain  number  of  grindings  as  the  number 
of  grindings  by  the  freezer  mill  increases,  which  could  be 
due  to  a possible  contamination  by  foreign  materials  during  the 
regrinding  process. 

By  judging  from  the  TEM  analysis  and  the  FIR  measurement, 
clustering  is  definitely  not  the  major  cause  of  the  anomalous  FIR 
absorption  and  the  sample  mixed  by  hand  which  contains  micron-size 
clusters  or  larger  did  show  an  enhanced  absorption.  Therefore,  the 
anomalous  FIR  absorption  is  indeed  the  small  particle  phenomenon. 
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Figure  3.7  Grinding  effect  on  FIR  absorption  of  Al-KCl  composites; 
the  arrow  indicates  the  result  of  Carr  et  al.'^ 
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Scaling  Behavior  in  e-i(o) 

As  a byproduct,  the  information  about  the  enhancement  of  the 

static  dielectric  constant  of  composites  can  be  analyzed  from 

interference  fringes  due  to  multiple  internal  reflections.  This  analysis 

is  most  easily  done  by  Fourier  transforming  corrected  (for  the  base  line, 
2 

K fv  ) the  absorption  coefficient.  This  transform  shows  a maximum 

located  at  x = 2nd  (n  = index  of  refraction  and  d = thickness  of  sample). 

An  example  calculation  is  shown  in  Figure  3.8.  The  curve  has  a 

Lorentzian  shape  because  the  interference  fringes  decays  exponentially 

as  the  frequency  increased  due  to  absorption.  The  volume  fraction 

dependent  behavior  of  e-|(o)  is  shown  in  Figure  3.9.  The  enhancement 

in  e-jCo)  with  increasing  small  particle  concentration  agrees  with  the 

42 

results  of  Grannan  et  al.,  described  in  Chapter  I. 


MIR  Data 

Absorption  Coefficient 

A dramatic  and  unexpected  change  in  the  absorption  coefficient 
was  observed  as  the  frequency  increased  above  the  FIR  frequency.  As 
shown  in  Figure  3.10,  there  is  a sudden  decrease  in  the  absorption 
coefficient  near  1,000  cm”\  after  which  the  absorption  increases 
quadratically  again  as  the  frequency  continues  to  increase.  However, 
the  MIR  absorption  strength,  is  much  smaller  than  in  the 

FIR  region.  The  MIR  data  were  fitted  to 


exp 


= ^expf'' 


(3.6) 
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Figure  3.8  An  example  of  obtaining  the  index  of  refraction  by  Fourier  transform  technique 
(a)  Base  line  corrected  absorption  coefficient  of  f = 0.016  A1  1-KCl  composite 
(d  = 0.055  cm)  and  (b)  Fourier  transformed  result  of  (a) 
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Figure  3.9  Scaling  behavior  of  Ci(0) 
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Figure  3.10  MIR  absorption  coefficients  of  A1-KC1  composites 
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The  measured  in  the  MIR  range  is  plotted  in  Figure  3.11  versus 
particle  size.  There  are  some  variations  in  K^^p  as  the  size  changes, 
but  it  is  quite  different  from  the  result  predicted  by  the  classical 
model  calculation. 

Note  that  the  classical  theory  given  in  Chapter  I explains  K^^p 
in  magnitude  if  we  use  a sufficiently  short  electron  mean  free  path, 

£ = VpT,  in  the  Crude  dielectric  function.  A fit  with  I = I K 
(t  = 4.95  X 10"^^  sec)  is  also  shown  in  Figure  3.11  as  the  solid  curve. 
The  satisfactory  results  with  an  unphysical  parameter  (1  ^ mean  free 
path)  are  rather  surprising  and  stimulates  the  speculation  about  the 
non-Drudian  possibilities  of  the  electron  behavior  in  a small  particle. 

The  MIR  absorption  data  of  the  Au  and  Ag  small  particles  are 
shown  in  Figure  3.12.  There  are  noticeable  differences  in  the  MIR 
behavior  of  Ag  and  Au  particles  from  that  of  aluminum:  The  MIR 

absorption  coefficients  of  Ag  and  Au  samples  do  not  seem  to  follow 
Equation  (3.5)  and,  in  addition,  there  are  many  narrow  structures  not 
seen  in  aluminum  small  particles.  Since  most  metal  oxides  have  a series 
of  optical  phonon  modes  in  the  frequency  range  between  300  cm'^  and 
1,000  cm  ^ and  the  small  particles  always  have  the  oxide  coating  on 
their  surface  which  was  introduced  for  the  electrical  isolation  between 
the  particles,  there  should  exist  the  optical  phonon  contribution  of 
the  oxide  in  the  absorption  mechanism  in  the  given  frequency  range. 

Oxide  Contribution 

The  effective  dielectric  function,  £p,  of  a small  particle 
constituted  of  a metallic  core  surrounded  by  an  oxide  coating  can  be 
written  as 
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Figure  3.11  Kgxp  in  MIR  frequency  vs.  particle  size;  the  classical 
calculation  with  t = 4.95  x 10"i^  sec  is  shown  as  the 
solid  line 
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Figure  3.12  MIR  absorption  coefficients  of  Au-KCl  and  Ag-KCl  composites 
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(3.6) 


~ ^ox  3e 


ox 


where  is  the  dielectric  function  of  the  oxide  coating  and 

? = ' - M'- 

Here  t is  the  oxide  thickness  and  a is  the  radius  of  the  particle 
(metal  plus  oxide). 


(resonance  frequency  and  oscillator  strength  for  each  mode)  were  chosen 


by  a factor  of  100  to  match  the  magnitude  of  the  broad  absorption.  Even 
so,  as  shown  in  Figure  3.13,  the  fit  to  the  experimental  observation  by 
using  the  parameters  and  Equation  (3.6)  gives  a location  for  the  peak 
in  the  calculated  plot  which  does  not  match  the  observed  peak.  Also,  an 
experimental  e for  evaporated  aluminum  oxide  (AlO  ) obtained  by 

OX  X 
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Eriksson  et  al.  was  used  for  comparison  and  the  result  is  shown  in 
Figure  3,13.  A fairly  good  agreement  (the  fit  is  about  30  cm"^  red 
shifted)  with  the  experimental  data  was  obtained,  which  supports  the 
electron  diffraction  result  that  the  oxide  coating  is  in  an  amorphous 
form  rather  than  a crystalline  form.  However,  as  shown  in  the  small 
insert  in  Figure  3.13,  the  high-resolution  experimental  data  shows  some 
weak  individual  phonon  peaks  of  aluminum  oxide  with  their  natural  width. 
This  observation  suggests  that  the  oxide  coating  consists  of  amorphous 
aluminum  oxide  with  small  fraction  of  crystalline  aluminum  oxide. 


In  the  case  of  the  aluminum  oxide,  the  parameters  for  e 


ox 


from  the  work  done  by  Barker. 
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But  the  linewidth  has  to  be  increased 


Another  possible  candidate  is  tungsten  oxide,  which  has  the 
higher  longitudinal  mode  than  the  aluminum  oxide.  However,  EDX  Analysis 
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A theory  tit  (dashed  curve;  to  MIR  data  by  using  parameters  given  in  Reference  b/ 
and  a data  fit  (dotted  curve)  with  experimental  value  obtained  in  Reference  68;  the 
small  insert  shows  a result  of  high  resolution  measurement 
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shows  no  direct  evidence  of  the  tungsten  contamination  in  AT  small 
particles  (see  Figure  2.5).  Also,  because  no  oxygen  gas  was  introduced 
during  the  aluminum  evaporation  process,  there  will  be  minimal 
oxidation  of  the  tungsten  filament  at  the  5.0  x 10  ^ Torr  background 
pressure  in  the  bell  jar. 

As  shown  in  Figure  3.12,  there  are  many  different  types  of  oxide 
contributions  to  the  absorption  coefficient  of  Ag  and  Au  composites. 

For  these  metals,  because  evaporation  took  place  in  approximately  10% 
oxygen  and  90%  of  Ar  inert  gas  mixture,  O2  also  oxidized  the  source 
material  (either  Mo  or  W).  Thus,  we  have  small  Ag  or  Au  particles  with 
an  oxide  coating  combined  with  WO^  or  MoO^-  A TEM  view  of  these 
samples  is  shown  with  the  corresponding  EDX  analysis  in  Figure  3.14. 
Therefore,  the  broad  absorption  in  the  frequency  range  between  300  cm~^ 
and  1,000  cm”^  is  mainly  due  to  optical  phonons  in  the  amorphous  oxide 
coating  and  the  absorption  characteristic  depends  on  the  type  of  metal. 

One  remark  on  Ag  and  Au  small  particle  data  is  that  because  now 
the  Ag  (Au)  smokes  are  the  mixture  of  Ag  (Au)  particles  with  oxide 
coating  and  the  source  material  oxide,  the  actual  metal  volume  fraction 
was  always  smaller  than  that  of  the  calculation  when  the  composite  was 
prepared.  This  point  will  be  clear  in  the  Maxwel 1 -Garnett  resonance 
data . 

Another  remark  is  the  nonquadratic  behavior  in  the  MIR  absorption 
coefficient  of  Ag  (Au)  particles  could  be  due  to  some  finite  (but  low) 
dc  conductivity  of  Ag  (Au)  oxide  as  will  be  discussed  in  Appendix  E. 
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Figure  3.14  A TEM  result  of  Au  small  particles  with  EDX  analysis; 
note  that  the  size  distribution  is  very  broad 
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Grinding  Effect-MIR 

The  evolution  of  the  MIR  absorption  coefficient  with  successive 
regrindings  of  the  random  composite  samples  is  shown  in  Figure  3.15. 

An  effect  of  grinding  on  the  MIR  absorption  coefficient  v;as  evident 
unlike  the  FIR  case.  The  MIR  absorption  data  shown  in  Figure  3.15(a) 
are  for  samples  mixed  by  hand.  Two  distinct  characteristics  can  be 
noticed:  One  is  that  the  frequency  dependence  on  the  MIR  absorption 

coefficient  is  not  quadratic  any  more  even  though  the  FIR  tail  still  is 
quadratic.  The  other  is  that  the  absorption  coefficient  is  not  linear 
in  the  volume  fraction  at  all.  The  nonquadratic  dependence  (rather 
square  root  dependence)  on  the  MIR  absorption  coefficients  can  be 
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explained  by  the  resistor-capacitor  model  worked  out  by  Sen  and  Tanner. 

In  Figure  3.15(b),  the  results  of  the  remixed  sample  (not  reground  from 

2 

the  pellet)  in  the  freezer-mill  are  shown.  The  evolution  from  v4Ij  to  oi 
dependence  on  the  MIR  absorption  coefficient  starting  from  the  lowest 
volume  fraction  can  be  clearly  seen.  However,  the  magnitude  of 
absorption  is  still  very  large,  but  decreased  from  the  samples  mixed  by 
hand.  Next,  the  plot  in  Figure  3.15(c)  shows  the  ordinary  absorption 
coefficient  after  regrinding  twice  from  the  samples  used  in  Figure 
3.15(b).  But,  the  volume  fraction  dependence  is  not  still  linear  and 
the  magnitude  of  the  absorption  is  decreased  again,  but  is  still  larger 
than  the  next  plot.  Figure  3.15(d).  The  results  of  the  final  stage  of 
the  grinding  procedure  is  shown  in  Figure  3.15(d);  the  absorption 
coefficient  now  shows  a linear  dependence  on  the  volume  fraction.  Also, 
the  same  effect  for  different  particle  sizes  is  shown  in  Figure  3.16. 


\ J 
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Figure  3,15  The  grinding  effect  on  MIR  absorption  data  of  A1  (A1  18) 
KCl  composites  (note  that  the  absorption  coefficients 
have  been  scaled  by 
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Figure  3.16  The  grinding  effect  on  MIR  absorption  coefficients  of 
A1  (a1  14)-KC1  composites 
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Therefore,  we  can  safely  conclude  that  if  the  absorption  coefficient 
is  linear  in  f,  then  the  sample  satisfies  the  random  homogeneity 
condition. 

Visible  Frequency  Data 

Now,  a couple  of  the  classical,  but  fundamental,  tests  can  be 
made  in  this  frequency  range.  For  aluminum  small  particles  interband 
absorption  can  be  studied;  for  Au  and  Ag  particles,  the  surface  plasmon 
mode  can  be  studied.  The  surface  plasmon  mode  is  not  seen  in  A1  because 
it  occurs  at  ca  = 6 X 10^  cm  ^ (7.6  eV). 

The  samples  were  diluted  10  times  from  f = 0.001  A1 , Ag, 
or  Au-KCl  composite  and  reground  4 times  for  the  optical  measurements. 

Interband  Transition  in  A1 

As  can  be  seen  in  Appendix  C,  the  parallel  band  absorption 

is  one  of  the  characteristics  of  simple  nontransition  metals.  Bulk 

aluminum  has  this  characteristic  absorption  at  12,000  cm  \ Optical 

measurements  on  the  diluted  aluminum  small  particles  composites  showed 

an  unexpected  result  that  the  location  of  the  interband  absorption  peak 

of  aluminum  small  particles  was  red-shifted  by  about  1,000  cm  ^ from 

that  of  a bulk  aluminum  as  shown  in  Figure  3.17. 

This  surprising  result  may  be  interpreted  as  a result  of 

expansion  of  the  lattice  constant  of  the  aluminum  small  particle. 

Because,  as  shown  in  Appendix  C,  the  absorption  is  caused  by  the 

pseudopotential  U2qq  of  aluminum  parallel  band  and  the  pseudopotential 

_2 

is  proportional  to  a where  a is  the  lattice  constant.  Therefore,  a 
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Figure  3.17  Interband  transition  of  aluminum  small  particles 
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change  in  the  lattice  constant  of  a few  percent  leads  to  the  red-shift 
in  the  interband  absorption  peak. 

In  addition,  there  is  another  broad  structure  developed  aroudn 


decrease  as  the  particle  size  gets  larger  (i.e.,  as  the  surface  to 
volume  ratio  decreases).  This  peculiar  structure  might  be  due  to  the 
broken  Bloch  symmetry  at  the  surface  which  has  been  ignored  in  the  bulk 
system  as  discussed  in  Chapter  I. 

Surface  Plasmon  Mode 

As  briefly  reviewed  in  Chapter  I,  another  characteristic  of  a 
small  particle  is  the  Maxwell -Garnett  resonance,  which  occurs  at  a 
frequency  related  to  the  plasma  frequency  by 


constant,  and  is  the  dielectric  constant  of  the  medium.  Here  the 

Crude  dielectric  function  was  used  for  e . 

m 

Because  Ag  and  Au  are  not  free  electron  metals,  the  screened 
plasma  frequency  enables  us  to  observe  the  surface  plasmon  when  Ag  or 
Au  particles  are  dispersed  in  a medium  with  dielectric  constant 
(KCl  was  used  in  this  study).  The  experimental  observation  of  the 
Maxwell -Garnett  resonance  of  Au-KCl  composites  is  shown  in  Figure  3.18 
and  the  resonance  of  Ag-KCl  composites  is  shown  in  Figure  3.19.  In 
order  to  fit  the  classical  theory  to  the  experiment,  the  volume  fraction 


20,000  cm  The  oscillator  strength  of  this  broad  peak  tends  to 


(3.7) 


where  cOp  is  the  bare  plasma  frequency,  is  the  core  dielectric 
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Figure  3.18  The  Maxwel 1 -Garnett  resonance  (surface  plasmon  mode)  of  Au-CKl  composite 
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3.19  The  Maxwel  1 -Garnett  resonance  of  Ag-KCi  composite 
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had  to  be  reduced  by  20%  which  confirms  that  there  were  less  small 

particles  than  estimated  by  the  weight  of  the  smoke. Also,  the 

observed  resonance  peaks  are  red  shifted  from  the  result  of  the  classical 

69 

theory.  This  problem  can  be  solved  by  oxide  coating  and  by 
remembering  that  there  was  extra  oxide  contribution  made  from  the  source 
material  and  this  source  metal  oxide  enhanced  leading  to  the 
red-shift. 

The  width  of  the  resonance  line  was  much  broader  than  the 
prediction.  Therefore,  the  order  of  10  A electron  mean  free  path  had 
to  be  introduced  to  explain  the  experiment.  As  shown  in  Table  1.1, 
many  similar  results  (broadening)  were  observed  by  other  workers. 
Therefore,  we  have  seen  a tendency  which  favors  the  nonmetal  lie 
possibility  in  the  small  particle  problem,  because  10  A meanfree  path 
is  unrealistic  if  we  assume  metallic  behavior. 

Nonmetal  lie  Possibility 

Now,  there  are  many  things  to  be  organized  at  this  point. 

Because  evidence  for  a behavior  in  small  particles  differing  from  the 
Drude  particle  point  of  view  have  been  observed,  some  tests  need  to  be 
done. 

Electric  or  Magnetic  Dipole? 

As  predicted  by  the  Maxwel 1 -Garnett  theory  in  the  low  frequency 
and  low  volume  fraction  limit  [see  Equation  (1.22)],  the  importance  of 
either  the  electric  dipole  term  or  the  magnetic  dipole  term  can  be 
determined  by  scaling  the  absorption  coefficient  by  CQ(a))^'^^  or  e^((o)^'^^. 
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respectively,  where  £q(w)  is  the  frequency  dependent  dielectric 

function  of  the  surrounding  medium.  Because  the  magnetic  dipole  term 

is  important  when  the  DC  conductivity  is  reasonably  high,  the 

conductivity  of  a small  particle  can  be  categorized  by  judging  from 

this  test.  The  results  are  shown  in  Figure  3.20.  Also,  an  example  of 
3/2 

scaling  by  is  shown  in  Figure  3.21  for  overall  frequencies. 

The  absorption  coefficients  scale  well  with  favoring  the 

electric  dipole  term  for  small  size  particles  not  the  magnetic  dipole. 

In  other  words,  the  intrinsic  conductivity  is  more  likely  poor 
according  to  this  analysis. 

Electron  Relaxation  Time 

The  fundamental  assumption  of  the  Drude  theory  is  that  there  is 
a constant  relaxation  time.  Now,  if  we  assume  that  the  dielectric 
function  of  a small  particle  is  the  Drude  dielectric  function, 

2 

Em(‘^)  = 1 - (3.8) 

0)  + iw/x 

and  calculate  the  x from  the  experimental  data  by  writing  the  absorption 
coefficient  such  as 

“exp^“^  " *^^^0’  “p’ 

(Here,  the  Maxwell -Garnett  effective  dielectric  function  with  only  the 
electric  dipole  term  is  used),  then  the  relaxation  time  x can  be 
calculated  by  inverting  the  function  such  as 
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Figure  3.20  FIR  absorption  coefficients  scaled  by  e and  e (co)  ' for  two  different 

particle  sizes 
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Figure  3.21  An  example  of  absorption  coefficient  before  and  after 
scaling  by  notice  that  the  optical  phonon 

frequency  range  of  KCl  is  replaced  by  the  data  taken 
from  A1 -paraffin  composites 
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T = G(e^,  f.  Wp,  CO,  ag^p) 


(3.10) 


where  the  Maxwell -Garnett  theory  was  used  to  describe  the  function. 
Then  the  particles  will  have  dc  conductivity. 


o 


ne^T 

m 


(3.12) 


This  conductivity  is  shown  in  Figure  3.22  for  <a>  = 235  K size  small 
particles.  What  we  expect  from  this  calculation  is  a constant  value 
when  we  plot  the  dc  conductivity  versus  frequency.  But,  the  results 
surprisingly  showed  that  the  behavior  of  the  small  particle  is  rather 
close  to  that  of  an  insulator!  Therefore,  we  have  the  definite  evidence 
of  the  non-Drudian  behavior  of  a small  particle. 
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Figure  3.22  of  an  aluminum  small  particle  calculated  from  the  experimental  data 


CHAPTER  IV 
THEORY 

Formalism 

In  this  chapter,  some  consequences  of  the  broken  translational 
symmetry  in  the  electron  motion  will  be  extensively  studied  within 
linear  response  theory  and  the  random  phase  approximation 
(self-consistent  field  theory).^^  72,49  theory  will  be 

compared  to  the  experimental  results  and  some  predictions  will  be 
made  for  the  experimental  observation. 

General  Theory:  Review 

When  a small  sphere  filled  with  electrons  is  placed  in  a 
semi-classical  electromagnetic  field,  the  response  of  electrons  to  the 
external  field  gives  a perturbation  to  the  density  operator  at 
equilibrium  such  as 

P = (4.1 ) 

to  linear  order  in  an  external  perturbation  H-j . By  choosing  a gauge 
in  which  the  scalar  potential  vanishes,  we  have 

^1  = ^ ("^-2) 

Here  t is  the  vector  potential  and  p is  the  electron  momentum.  The 
unperturbed  density  operator,  p^,  satisfies 
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P°  |i>  = f(E.)  |i>  E f.  li>  (4.3) 

where  f^.  is  the  Fermi-Dirac  distribution  function  and  | i>  are  the 
eigenstates  of  the  unperturbed  free  particle  Hamiltonian  Hq  with 
eigenvalue  E^. . Then  the  Liouville  equation 

iTi  If  = [H,p],  H = Hq  + (4.4) 

can  be  linearized  by  ignoring  the  second-order  term  in  From 
Equations  (4.1)  and  (4.4),  we  have 

iTi^<irp^  Ij>  = <il  [H^,p°]  lj>  + <i|  [Hq,p^]  |j>  (4.5) 

where  the  higher  order  term  [H-j,p^]  has  been  ignored  and  the 
unperturbed  contribution 

i?i  ^<i|  p°  |j>  = <i[  [Hq,p°]  |j>  (4.6) 

has  been  subtracted.  With 

(4.7) 

making  the  ansatz  that  <i | p [j>  has  the  same  time  dependence  as 

^(r,t),  and  using  the  in  prescription  to  switch  adiabatical ly  on  the 

perturbation  from  t = -<»  (throughout  this  chapter,  lim  is  implicit), 

rr^O 


we  have 
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<i| 

By  using  Equations 
<i| 


' E.  - E.  - Tto  - in 

' vi 


(4.8) 


(4.2)  and  (4.7),  we  can  rewrite  Equation  (4.8)  as 


IJ>  = S ■ in  IJ> 

^ 3 


,M<i|  e’5-^  |j>] 


(4.9) 


with  E.j  = E,  - Ej. 


The  induced  current  density  (r,t)  may  be  derived  from 


T rO“^(^)x  ^'t(^)l 

Jind  ('"’■t)  = Tr  {p  j'  ' + p o' 


(4.10) 


where 

= (4.11a) 

= ^«(7-7op) 


and  r^p  are  the  momentum  and  position  operators,  respectively.  By 
using  the  relation  for  an  operator  in  a statistical  ensemble 


Tr(pJ^)  = I <i\  p |j>  <j1  |i> 

ij 


(4.12) 
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we  have 


'^ind 


I fi  <'l  - *op> 

TJ  ^ 


- I <■'1  |j><j|  h‘>  (4.13) 

ij 

Carrying  out  the  Fourier  transformation  on  Equation  (4.13),  and  using 
Equation  (4.9),  we  obtain 


(l^»co) 


'ind  I <’■! 


i (q-t) *r 


Ii> 


2 f.  - f. 

_§_  y 1 J 
,,  2 .4  E . . - Tioj  - in 

Vm  c ij  ij 


<i|  e’'’”'  J(5,<o)  • f \i> 


+ 1 Tiq  • !^(q,a))  <i  | |j> 


<3- 1 


n> 

(4.14) 


where  V is  the  volume  of  the  system  and  t is  the  Fourier  component 


vector. 
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By  assuming  homogeneity  in  the  system,  i.e.,  that  there  is  a 

response  only  at  the  wave  vector  of  the  external  field  (i.e,,  that 

q = t),  and  by  taking  the  q ->■  0 limit)  dipole  approximation),  as 
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appropriate  for  a small  particle,  we  find 


■^i  nd 


(O)) 


~Y  A(w) 


2 E . . 
Vm  c ij  ij 


- f. 

J 

Tio)  - iq 


• <il  p |j><j|  p 1 i>] 


(4.15) 


where  ^ f^. , the  total  number  of  electrons  in  the  system. 

Therefore,  we  can  derive  the  conductivity,  a(co),  and  the  dielectric 
function,  e(o)),  by  using  the  following  formulae  with  the  Coulomb  gauge. 

t((o)=^A(o))  (4.16) 


^ind 


(m)  = a((o)  t(co) 


(4.17) 


and 


;((,)  = 1 + 4TmM 

u 


(4.18) 


:(<») 


From  Equations  (4.15)  and  (4.17),  we  have 


a(w)  = i 


^2  ^2 


f.  - f. 

1 JL 


V ' _2  A E..  - tiw  - in  P 
Vm  (0  ij  ij 


(4.19) 


Here  the  vector  notation  in  p represents  the  direction  of  the  momentum 
along  the  applied  field.  Also,  the  dielectric  function  e(o))  can  be 
found  from  Equations  (4.18)  and  (4.19): 


2 


“r. 

e(w)  = , 1 - "2 


K^l  p IJ>I^ 


CO 


(4.20) 


behavior  term  in  the  square  brackets  and  a correction  term  due  to  the 
broken  symmetry  in  the  third  term  of  Equation  (4.20).  This  correction 
term  will  vanish  in  the  a ->  «>  limit.  This  limit  describes  the  free 
electron  system  without  relaxation  (i.e.,  t t = electron 
relaxation  time. 

It  is  normally  assumed  that  n goes  to  zero  after  the 
integration.  This  is  a prescription  included  in  the  RPA  to  prevent 

49 

the  divergence  of  the  time  integral.  In  contrast.  Wood  and  Ashcroft 
assumed  an  actual  physical  meaning  to  the  in  term.  They  identified 
n with  ti/x,  T being  the  usual  (Drude-model ) relaxation  time. 

A different  interpretation,  however,  will  be  made  in  this 
dissertation,  which  uses  the  identity 


(4.21) 


where  P stands  for  the  principal  part.  Then  we  can  explicitly 
separate  the  real  part  and  the  imaginary  part  of  the  dielectric 
function  such  as 
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e-|  (co)  = Re  e(o))  = 1 


4TTe^  V 

2 " mJ  2 A 
CO  Vm  CO  ij 


<i|  P |j>  ^ 


(4.22) 


and 


£2^“) 


= Im  e(co)  = I (f.  -f.)I<il  p |j>1^6(E..  -Tico) 

Vm  CO  ij 


(4.23) 


Also,  from  Equation  (4.18),  we  can  find  the  conductivity  such  as 


cj(co)  = -^  [1  - e(co)]  = o^(co)  + io2(co) 


(4.24) 


Equations  (4.22)  and  (4.23)  describe  divergences  in  e(co)  at 

Tico  = E..  = E.  - E.  and  give  infinitely  sharp  absorption  lines  in  e^. 

In  practice,  we  never  see  such  sharp  absorption  lines  because  of 

broadening  effects.  For  a small  particle,  a quantum  mechanical 

broadening  needs  to  be  considered,  as  discussed  in  Chapter  I and 

Appendices  A and  B,  since  the  quantum  mechanical  role  of  the  surface  is 

not  that  of  a scattering  center,  but  in  determining  the  eigenstates  of 

bound  electrons.  Because  the  random  irregularity  of  the  surface  plays 

an  important  role  (level  correlation),  averaging  a physical  quantity 

74 

over  impurities  and  imperfections  on  the  surface  is  the  most 
important  contribution  to  the  broadening. 
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Small  Particle  and  Matrix  Element 
The  crucial  point  in  the  RPA  small  particle  theory  is  the 
evaluation  of  the  matrix  element 

= l<i|  P |j>l^  (4.25) 

The  momentum  matrix  elements  given  by  Equation  (4.25)  can  be 
transformed  into  a position  matrix  representation,  which  is  more 
appealing  in  the  case  of  a small  particle  because  the  effect  of  broken 
symmetry  can  be  directly  visualized,  by  making  use  of  the  Ehrenfest 
theorem 

^ <SX>  = - 1 [f2,H]  (4.26) 

Thus,  we  have 

<il  p 1j>  = -i  ^<i|[r,HQ]  1j>  = -i  ^ E^.^.  <i  | r 1j>  (4.27) 

Here  the  vector  notation  represents  the  direction  along  the  applied 
field  again.  Therefore,  the  momentum  matrix  elements  will  be 
replaced  with  the  corresponding  position  matrix  elements  whenever  we 
need  a clear  physical  description  of  e(o))  and  a(co). 

The  matrix  elements  given  by  Equation  (4.27)  can  be  evaluated, 
in  principle,  if  we  know  the  eigenstates  and  eigenvalues  of  a given 
Hamiltonian  from  the  solution  to  the  Schrbdinger  equation.  However, 
finding  the  solutions  to  the  Schrbdinger  equation  is  impossible  in  our 
model  because  of  the  random  nature  in  the  boundary  condition.  Hence, 
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a reasonable  approximation  needs  to  be  made  to  evaluate  the  matrix 
elements. 

Low  Frequency  Limit 

Now,  because  the  matrix  elements  generally  depend  on  the  energy 
eigenvalues  of  given  eigenstates,  for  frequencies  below  a minimum 
energy  gap  of  the  system,  we  can  approximate  the  matrix  elements  by 
considering  electrons  on  the  Fermi  surface  such  as 

|<i|  r |j>|^  E \r..\^  = 1 Ir.j.l^  (4.28a) 

and 


^a^ 
5 ® 


(4.28b) 


y 

Here  we  have  taken  the  spatial  average  on  k.jjl  • 

For  frequencies  above  the  minimum  energy  gap,  we  expect  that 

electrons  in  a small  particle  will  recover  some  "free"  electron  nature. 

Thus,  we  can  treat  the  problem  semiclassical ly  by  employing  the 
75 

Shapoval  method  when  the  De  Broglie  wave  length  of  electrons  do  not 
exceed  the  roughness  of  the  surface  (i.e.,  the  diffusive  scattering 
case) . 


Starting  from  the  Heisenberg  picture,  the  average  of  two 
physical  observables  f(t)  and  g(t  + t')  in  a state  | i>  at  a time 
difference  t'  can  be  calculated  as 


104 


<f(t)g(t  + t')>.  = <iff{t)g(t  + t')li> 


I <i|f(t)|jXjlg(t  + f )|i>  (4.29) 


By  changing  Equation  (4.29)  to  the  Schrbdinger  picture,  we  have 


-ico.  .t' 

<f(t)g(t  + t')>^.  = I e <ilf| jXjjgfi> 

where  we  have  defined  co . . as  w - • h For  f = g = r,  we  have 

I J • J II 


(4.30) 


-iw..f  p 

<r(t)r(t  + t')>^.  = ^ e ^ ' 

j 


(4.31) 


Now,  by  changing  the  summation  to  an  integration,  we  can  rewrite 
Equation  (4.31 ) as 


<r(t)r(t  + t'  )>^.  = Tiv(Ep) 


e-'"*'  do, 


where  v(Ep)  is  the  density  of  states  at  the  Fermi  level 


(4.32) 

Vmpp 


^^^F^  " 2,.3’ 


Pf 


= Fermi  momentum 


TT  t)' 


. After  the  Fourier  transformation,  we  get 


’■ij< 


2ntlv(Ep)  , 


<r(t)r(t  + t' )>^  e 


ioii-f 


dt'  (4.33) 


Next,  the  quantity  in  the  angular  brackets  of  Equation  (4.33)  can  be 
averaged  by  taking  into  account  the  classical  trajectory  of  electrons 
in  the  diffusive  scattering  limit  as 
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<r(t)r(t  + t'  )>.  = -2Y 


T 

-T 


dt  r(t)r(t  + t' ) 


(4.34) 


The  autocorrelation  function,  <r(t)r(t  + t')>  vanishes  when  t' 
coincides  with  a/vp,  which  is  the  minimum  time  interval  to  have  a 
surface  scattering.  Also,  since  the  autocorrelation  function  is  an 
even  function,  we  can  rewrite  Equation  (4.33)  in  the  form  of  a cosine 
Fourier  transform; 


TTTlv(Ep) 


<r(t)r(t  + t' )>. 


0 


cos(w.  .t' 

* J 


)dt' 


(4.35) 


For  very  low  frequencies  (co  « Vp/a),  we  can  ignore  the  frequency 
dependence  of  the  matrix  elements  because  cos(o)t')  is  close  to  unity. 
Thus,  the  matrix  elements  are  independent  of  the  energy  states  in  this 
limit.  Therefore,  we  can  set  the  matrix  elements  to 


2 1 


ij  'o  TrTiv(Ep)  _ 


<r(t)r(t  + t'  )>^.  dt 


1 1 


TTTiv(Ep)  2T 


dt 


dt'  r(t)r(t  + t' ) (4.36) 


-T  0 


The  integration  given  by  Equation  (4.36)  can  be  done'  to  give 


2 139  a' 


i j ' 0 150  27rTiv(Ep) 


(4.37) 
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High  Frequency  Limit 

In  the  limit  of  high  frequencies  (co  » Vp/a),  a quasi-continuous 

behavior  of  electrons  in  a small  particle  is  anticipated.  In  other 

words,  unlike  the  low  frequency  limit,  the  details  of  discrete  level 

structures  become  unimportant.  However,  in  our  small  particle  model , 

surface  roughness  is  expected  still  to  play  an  important  role.  This 

point  will  be  clear  later  in  the  surface  plasmon  broadening  section. 

For  a sphere  with  a smooth  surface,  we  can  easily  find 

eigenstates  of  the  Schrbdinger  equation  in  terms  of  spherical  Bessel 

functions  (see  Appendix  A and  Reference  48).  But,  working  with  these 

eigenstates  to  evaluate  the  position  matrix  elements  is  cumbersome. 

Furthermore,  it  is  impossible  to  determine  the  eigenstates  for  a 

48 

sphere  with  a rough  surface.  However,  Kawabata  and  Kubo  simplified 
the  situation  for  a smooth  sphere  by  introducing  a surface  potential, 
U(f),  as  follows: 


<i|  r I j>  = <i|  [H,r]  |j> 

ij 


<i 

m 

ij 


5U(r) 


->• 

8r 


I j> 


(4.38) 
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from  the  commutation  relations 


[H,r]  = P and  [H,p]  = in  (4.39) 

dr 

This  idea  can  be  used  in  our  small  particle  model  by  treating  U(r) 
as  a position  dependent  random  surface  potential  to  take  into  account 
the  effect  of  surface  roughness.  However,  further  systematic 
investigation  needs  to  be  done  in  the  future. 


Calculation  of  e((o) 

From  Equations  (4.22),  (4.23),  (4.24),  and  (4.27),  we  have  the 
desired  expression  for  e(co)  as  follows: 


4iTai  (co) 

e(co)  = £■]  (co)  + 


(4.40) 


where 


e-|(co)  = 1 - 
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CO 
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ij 


fi 
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- Itoo 


'ij 


hco 


|2 


(4.41) 


and 


(co)  = 


Tre  CO 
V 


1 

ij 


ij 


(E,j  - fw) 


(4.42) 


e(co)  given  by  Equation  (4.41)  can  be  rearranged  in  two  equivalent 
forms  to  give  a better  convergence  in  a series  expansion  or  a clearer 
physical  meaning  by  interchanging  indices  i and  j in  the  summation 
in  Equation  (4.41 ) to 
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or 
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£l(w)  = 1 - 
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(4.43) 
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Level  Correlation 

The  most  interesting  physical  quantities  in  a small  particle 
are  correlation  functions.  The  energy  correlation  function  of  a small 
particle  R(oj),  defined  as 

R(oj)  = < ! ^ y (f,  - f,)  6(E,,  - T«o)>  (4.45) 

4rk»[v(Ep)]^  IJ  J 


has  a direct  connection  to  the  response  to  an  external  electromagnetic 
field,  because  it  gives  the  probability  of  finding  two  energy  levels 
a distance  fiw  apart.  Here  the  angular  brackets  in  Equation  (4.45) 
represent  the  average  over  all  impurities  and  imprefections  in  the 
system.  This  averaging  process  is  necessary  to  include  the  quantum 
broadening  as  discussed  earlier  in  this  chapter. 
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Generally,  the  evaluation  of  averaged  physical  quantities  like 
Equation  (4.45)  is  difficult  because  it  is  impossible  to  find  exact 
solutions  of  the  Schrbdinger  equation  in  such  systems.  Perturbation 
techniques  can  be  used  in  limited  cases,  when  the  Hamiltonian  can  be 
expended  in  a series  over  sufficiently  weak  random  potentials. 

Another  method  developed  to  calculate  such  averaged  quantities 
is  the  replica  method, applied  first  to  the  spin  glass  problem, 
which  replaces  the  random  system  by  n subsystems  identical  to  the 
original  and  uses  the  identity 


fi,n  X = lim 
n->0 


(4.46) 


Thus,  one  can  carry  out  the  averaging  procedure  before  taking  the 
n 0 1 imit  (e.g. , for  a free  energy,  F = - kgT  Jin  Z,  we  can  calculate 
the  average  value  of  Z first  and  then  take  n 0 limit).  However, 
there  are  many  difficulties  in  applying  this  method  because  analytical 
continuation  or  an  integer  condition  for  n in  the  standard 

approximation  usually  fails  in  the  n 0 limit. 

78 

Recently,  Efetov  evaluated  the  two-level  correlation  function 
given  by  Equation  (4.45)  which  can  be  rewritten  in  terms  of  the 
retarded  (-)  and  advanced  (+)  Green's  function,  G",  such  as 


R(o)) 


8it^w[v( 


E IT 


<Re 


dE 


I 


[f(E+Tio))  -f(E)]G|(E) 


{Gt(E+tico)  - G;j(E+Tia))»  (4.47) 

J J 
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by  applying  the  supersymmetry  method. The  main  advantage  of 

using  this  method  is  that  the  Green's  function  can  be  reexpressed  in 

the  form  of  functional  integrals  by  employing  the  supervector  and 

the  supermatrix  algebra  and  the  new  Green's  function  does  not  have  the 

weight  denominator  which  appears  in  the  ordinary  form. 

Now,  we  can  take  the  average  on  Equation  (4.47)  at  the 

beginning  on  the  product  of  new  Green's  functions  and  this  procedure 

78 

requires  a fair  amount  of  mathematical  art.  Surprisingly,  Efetov 
found  results  which  exactly  coincide  with  those  derived  on  the  basis 
of  Dyson's  phenomenological  theory  for  three  different  ensembles, 
shown  in  Appendix  B,  by  using  a zero-dimensional  supermatrix  a-model . 

With  the  aid  of  this  important  consequence  of  the  new 
calculation  of  the  two-level  correlation  function,  we  can  evaluate 
averaged  quantities  over  surface  imperfections  and  impurities. 

Low  Frequency  Limit 

In  the  low  frequency  limit  where  the  discreteness  of  energy 

levels  and  their  correlation  is  important,  a subtle  distinction  is 

necessary  to  obtain  the  correct  low  frequency  form  because  of  the 

2 2 

diamagnetic  current  contribution  which  results  in  Wp/co  term  in 
Equation  (4.22). 

If  we  can  choose  a minimum  gap  energy,  one  may  expand 

the  formula  given  by  Equation  (4.44)  for  the  frequency  range  below 
AEmin  to  obtain  the  low  frequency  form  such  as 
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(4.48) 


where 


E 1 + 


Sfre 


1 


(4.49) 


We  can  show  that  the  optical  effective  mass  (m*)  defined  by  Equation 
(4.49)  becomes  infinite  (i.e.,  there  is  no  continuum  contribution)  in 
a truly  discrete  quantum  system  for  frequencies  smaller  than  the  gap 
by  applying  the  Thomas-Reiche-Kuhn  sum  rule 


(4.50) 


and  ^ f.j  = Ng.  Thus,  the  optical  effective  mass  ratio  vanishes  in  a 
small  particle  for  the  same  reason.  Therefore,  we  may  rewrite  the 
general  result  given  by  Equation  (4.41)  as 


Cl(co)  = 1 - 
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Ir. 


> 


(4.51) 
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for  frequencies  below  the  gap  by  following  the  argument  given  in 
Equation  (4.48).  Here  the  angular  brackets  represent  the  average  over 
all  impurities  and  imperfections  on  the  surface.  Also,  the  same 
averaging  process  needs  to  be  considered  for  ai(u)  to  give 


a^(co)  = 


^2 

Tie  CO 


<I 

ij 


Ir: 


-/6{E..  -r^)> 


(4.52) 


At  this  point,  we  have  confirmed  the  Kubo  formula  derived  from  the 
fluctuation-dissipation  theorem  which  leads  to 


a-|  (co)  = Re  a(co) 


dt 

0 


3 

dX  e”''“^  4(0)  j(t  + iX)> 

0 


(4.53) 


Here  6 = 1/kgT  and  j^(t)  represents  the  current  operator  in  the 

Heisenberg  picture.  Equation  (4.52)  can  be  recovered  immediately  from 

81 

Equation  (4.53)  after  some  algebra. 

Now,  the  static  part  of  the  electric  susceptibility  defined 
as 


Xp  = lim  xM  = lim 
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4ir 


(4.54) 
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(4.55) 
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can  be  evaluated  as 


x2  = 


3 <'^p) 

f ^ 
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00 
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oj 

(4.56a) 


or 


5tt^ 


(4.56b) 


where  a^=  Bohr  radius  = — 2»  “ m v^.,  and  y is  the  specific  heat  mass 

ratio  (m  /m).  [N.B.,  y is  different  from  the  mass  ratio  defined  by 

Equation  (4.49).]  We  have  used  the  following  relations  to  obtain 
Equation  (4.56)  from  Equation  (4.55): 


fi  - f-  -6(Ep  - E.)  E.j 


v(Ep)  = 2 I 6(Ep  - E.) 

j 


(4.57a) 

(4.57b) 


and  Equation  (4.28). 

By  taking  into  account  the  statistical  nature  of  the  energy 
gap  of  a small  particle  as  discussed  in  Appendix  B,  we  obtain  the 
frequency  dependent  susceptibility, 


114 


XgU) 


<I 


ej 


ij  E..  - (Hco) 


AE 


2l 


^v(Ep) 


v{Ep) 


dE 


P(E) 


E^  - (Tico)' 


dx 


P(x) 


0 x' 


flO) 


AE 


2 


X = E/AE 


where  P(E--)  is  the  level  distribution  function  which  depends  on  c 

* vJ 

statistical  ensemble  in  the  system.  These  P(E^j)  have  been  shown 
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(symplectic  ensemble) 


For  the  orthogonal  ensemble,  for  example,  we  find 
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Xg  (w)  = X 

e 


1 - ^ 


9 IT 


AE 


AE 


Ei 


W 

AE 


(4.58) 
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(4.62) 


(4.63) 
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with  the  aid  of  the  formula 
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+ I (k  - 1)!  (-3)"'''  y"''  (4.64) 

k=l 
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from  Gradshteyn  and  Ryzhik.  Here  Ei(x)  is  the  exponential  integral 
defined  as 


Ei(x) 


(4.65) 


But,  one  difficulty  arising  in  the  analytic  evaluation  of  Xg(w)  is  the 
limitation  on  the  range  of  x in  Ei(x)  because  the  level  distribution 
function  covers  the  frequencies  up  to  Tico/AE  25  for  a smooth 
connection  to  Xg(^i>)  for  Tico  > AE.  Therefore,  Xg(^*j)  was  evaluated 
numerically  in  this  work. 

Also,  the  real  part  of  the  conductivity,  a^(oj),  given  by 
Equation  (4.52)  can  be  calculated  as  the  same  way  done  in  the  Xg(w) 
calculation  to  give 
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Therefore,  for  Tico  < AE,  we  have 
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Here 

= |-  and  y-j  = ^ orthogonal  ensemble  (N  = 1) 

Co  = ^ and  yo  = - for  unitary  ensemble  (N  = 2) 
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for  symplectic  ensemble  (N  = 4) 


One  interesting  prediction  of  this  calculation  is  that  the  dc 
conductivity  defined  as 


0dc  = 1 ini  o,  (w)  (4.69) 

or>0 

vanishes  for  a small  particle.  Classically,  this  point  can  be  viewed 

as  a finite  q case  (note  that  the  Drude  dielectric  function  is  a 

limiting  case  of  q = 0).  Therefore,  a small  particle  is  an  insulator!* 

For  frequencies  above  the  gap,  we  cannot  use  the  low  frequency 

form  of  e-]  (oj)  given  by  Equation  (4.51)  directly  because  the  role  of 

the  diamagnetic  current  becomes  important  in  this  frequency  range. 

2 2 

Thus,  by  rewriting  the  cOp/co  term  in  Equation  (4.43)  as 

*The  term  insulator  means  the  vanishing  dc  conductivity  at  zero 
temperature. 
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(4.70) 


we  obtain 


(4.71) 


without  loss  of  general ity.  Note  that  Equation  (4.71)  is  identical  with 
Equation  (4.51)  except  that  the  valid  frequency  range  is  above  the  gap 
and  the  position  matrix  elements  are  different  in  this  frequency 
range.  In  other  words,  the  contribution  of  the  static  part  of  the 
susceptibility  should  vanish  for  this  frequency  range.  This  result 
provides  a correct  physical  description  of  a small  particle  because 
there  will  be  no  contributions  of  transitions  across  the  gap  to  the 
susceptibility  at  frequencies  above  the  gap  like  ordinary  nonmetallic 
systems. 

Because  the  electrons  in  a small  particle  are  expected  to 
recover  some  free  electron  character  for  frequencies  above  the  gap,  we 
have 


(4.73) 
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from  the  definition  of  R(o))  given  by  Equation  (4.43)  and  the  matrix 
elements  given  by  Equation  (4.37). 

For  £i(co)  above  the  gap,  we  obtain 
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= 1 + 4-rr 
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(4.75) 


because  the  effects  of  the  averaging  process  are  now  absorbed  into  the 
level  correlation.  The  dimensionless  function  F(w)  defined  as 


F(u)  = (fio))^  v(Ep) 


oo 


(4.76) 


is  identical  to  the  real  part  of  the  function  A(co)  defined  by  Gor’kov 

7 85 

and  Eliashberg'  [however,  note  that  Devaty  and  Sievers  reevaluated 

the  function  A(co)  and  corrected  errors  in  the  original  calculation]. 

F(w)  for  three  different  ensembles  are  given  by 


F'"'t(x)  = 
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(4.77) 
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and 


psympl(^)  , 2 . sin_^  , S,(y)  ^ (cop- 


(4.79) 


where 
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cos  t I r*/  \ sin  t II 
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A general  behavior  of  e-|(oj)  and  a-|(a))  for  three  different 
ensembles  is  shown  in  Figure  4.1.  One  remarkable  result  of  this  theory 


Gor’kov  and  Eliashberg  case  which  will  be  shown  in  the  next  section. 

Note  that  the  static  part  of  the  dielectric  function  is  effectively 

canceled  out  as  the  frequency  increases  across  the  gap  frequency  and 

the  dielectric  function  converges  toward  1 and  makes  a smooth 

connection  to  the  result  of  Equation  (4.75)  at  frequencies  away  from 

the  gap.  Also,  note  that  the  behavior  of  £i(cd)  is  quite  similar  to  that 
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of  an  ideal  insulator  at  the  visible  and  ultraviolet  frequency  ranges. 

In  other  words,  the  same  physical  phenomena  occuring  at  high  frequencies 
in  a nonmetal  lie  system  have  been  brought  down  to  FIR  or  microwave 
frequencies. 

The  particle  size  dependence  of  e.|(a))  and  the  corresponding 

a.|  (o))  are  shown  in  Figures  4.2  and  4.3,  respectively.  Notice  that 

2 

the  small  particle  conductivity  follows  a co  dependence  above  the  gap 

frequency  and  a small  particle  has  intrinsically  poor  ac  conductivity 

59 

which  resembles  exactly  the  three-dimensional  Anderson  insulator! 


is  that  £.|(a))  naturally  converges  to  (1  + 4itx°)  as  co  0 unlike  the 
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Figure  4.1  General  frequency  dependence  of  C](a))  and  O]  (co)  near  the  gap  for  a 20  A alumi 
small  particle 
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Additionally,  a comparison  among  a = 20  A A1 , Au,  and  Sn  small 
particles  for  three  different  ensembles  is  made  in  Figures  4.4  and 
4.5  for  e-|(o3)  and  a-](co),  respectively. 

Gor’kov  and  Eliashberg  Theory 

Gor’kov  and  Eliashberg  (GE)^^  applied  the  idea  of  energy  level 
correlation  in  a random  system  to  the  small  particle  problem  and 
obtained  the  electric  susceptibility,  Xg’  evaluating  the 
expectation  value  of  the  electric  dipole  moment  operator,  ^ = er, 
directly  as 


<3>  = <ef>  = V Xg  <&>  (4.80) 

with  a perturbation  H-j  = By  using  the  same  approximation  used 

in  the  general  theory  review  section  of  this  chapter,  we  find  the 
electric  susceptibility  as 

Xg  \t\^  = <^>  • t/\l  (4.81) 

\t\^  (4.82) 

Therefore,  from 
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Figure  4,4  e-,(v)  for  different  types  of  metal  small  particles  for  three  different  ensembles 
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we  find 


e(w)  = 1 


4ire 
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ij 
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(4.84) 


which  is  essentially  the  same  as  Equation  (4.51).  At  this  point,  GE 
applied  the  two-level  correlation  function,  R(co),  conjectured  by 
Dyson^  to  change  the  summation  to  an  integration  after  separating  e(co) 
into  static  and  dynamic  parts  and  found 


(4.85) 


where 


and 
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Here  the  function  A(m)  is  defined  as 


A(o))  = (Tiw)^  v(Ep) 


■ R(E)  dE 
E^  - (Tiw  + in)^ 

— oo 


(4.86) 


Note  that  the  real  part  of  A(w)  is  identical  to  the  dimensionless 
function  F(w)  introduced  in  Equation  (4.75).  The  imaginary  part  of 
A(w)  can  be  found  as 
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and 


A°>-‘{x)  = 2x  (2x) 

= (2x)  R°''’^(x) 

^unit(x)  , 2^  . 

= (2x)  R^"'^(x) 


- Si(x) 


nl 

2]  dx 


sin  x' 


(4.86a) 


(4.86b) 

(4.87a) 

(4.87b) 


A2^'"'^'(y)  = y 


y R 


sympl 


(y) 


where  y = 2x  = — and  A(aj)  = A-|(co)  + iA2(o)) 
Therefore,  we  have 


a-|(w)  = Q A2(o)) 


(4.88a) 

(4.88b) 


(4.89) 


and 

£l(a))  = 1 + 4ttXq^  + 4ttQ  Ai(co)  (4.90) 

GE  results  given  by  Equations  (4.89)  and  (4.90)  are  shown  in  Figure 
4.6.  Note  that  Equation  (4.90)  does  not  have  the  correct  physical 
description  in  a sense  that  the  static  part  of  the  electric 
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Figure  4.6  £](w)  and  a-jCcj)  of  Gor’kov  and  Eliashberg  theory;  note  that  the  small  insert  in  e-j(a)) 
plot  shows  the  oscillatory  nature  arising  from  the  two  level  correlation 
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remains  unchanged  as  frequency  increases  across  the  gap.  Also, 
notice,  therefore,  that  0-1(00)  does  not  have  the  correct  behavior  across 
the  gap.  All  these  discrepancies  are  due  to  the  lack  of  information 
about  the  diamagnetic  current  contribution  in  the  induced  dipole 
moment  approximation.  However,  for  the  frequencies  far  below  the 
gap,  the  GE  theory  qualitatively  explains  the  electromagnetic 
properties  of  a small  particle  because  the  low  frequency  tail  of  the 
level  correlation  function  [R(E)j  has  the  same  behavior  as  the  level 
distribution  function  [P(E)]. 

High  Frequency  Limits 

In  the  limit  of  high  frequencies,  we  expecr  that  the  level 
correlation  loses  its  meaning  and  the  system  can  be  treated  in 
a quasi-continuum  limit.  However,  the  calculation  of  the  position 
matrix  elements  with  random  surface  potential  (random  superposition 
of  individual  atomic  potential)  is  important  because  the  origin  of 
broadening  has  to  come  from  the  random  nature  of  the  surface.  But, 
at  the  present  time,  no  systematic  method  has  been  developed  and 
research  needs  to  be  done  on  this  problem  in  the  future. 

The  Classical  (Bulk)  Limit 

The  a •>  °°  Limit 

As  the  geometrical  size  of  a particle  increases,  the  average 
energy  gap  will  approach  zero  and  the  cutoff  frequency  co^  = Vp/a  in 
the  low  frequency  matrix  elements  will  also  become  zero. 

For  the  conductivity  a-|(a)).  Equation  (4.68)  will  become  a 
delta  function  as  AE  ^ 0 from 
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6(x)  = lim  -3  e ^ 


(4.91) 


This  is  the  case  of  metal  without  relaxation  as  given  by  Equation  (1.1). 
Also,  £.|(a))  given  by  Equation  (4.44)  converges  to 


as  it  should  be  in  the  bulk  limit  of  electron  gas  without  relaxation. 
The  two  limiting  case  of  a.j((o)  and  e.|(a))  in  our  calculation  is 
remarkable  because  Equations  (4.68)  and  (4.44)  show  the  correct 
behavior  unlike  the  GE  calculation. 

Now,  if  there  exists  a finite  electron  mean  free  path  (£),  a 
small  particle  is  expected  to  recover  the  classical  (Drude) 
conductivity  from  the  small  particle  conductivity  as  the  size  of  a 
particle  exceeds  the  electron  mean  free  path.  To  make  this  argument 
more  transparent,  the  position  matrix  elements  need  to  be  reexamined 
in  this  limit.  Because  there  exist  random  scattering  centers  for 
electrons  (or  electrons  have  a finite  lifetime,  t = il/vp),  an 
exponentially  decaying  behavior  in  the  correlation  function  may  be 
assumed.  Thus,  the  electron  position  correlation  function  given  by 
Equation  (4.31)  can  be  rewritten  as 


2 


(4.92) 


<r(t)  r(t  + t')>^.  = I |<i  I r(0)  |£>|  ^ e~^^ 


(4.93) 
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Hence,  the  position  matrix  element  can  be  calculated  from  Equation 
(4.33)  to  give 


TrtivlEp)  1 


— y 

T\/,\  L 


riu  i 


(4.94a) 


1 T 

03[1  + (<0j,T)^] 

' u 


(4.94b) 


where  \ (E.  - E^) ]<i  | r(0)  | O]  ^ = - ^ and  E^.  - E^  'v  Tico  have  been  used. 

Therefore,  the  Drude  conductivity  can  be  recovered  by 
rearranging  Equation  (4.42)  in  terms  of  the  absorption  and  emission  of 
radiation  parts  such  as 


ai(o))  = 


^2 

ire  03 


I 

Tj 


Ir. 


[6(E, 


+ no))  - fi(E^-j 


- no3)] 


(4.95) 


and  by  changing  the  summation  about  j to  integral  as 


I ^ -hv(Ep) 

j 


00 


do) 


—00 


(4.96 
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we  have 


’ 2mV 


do) 


f \ 

CO  . . 

J1 


1 + (co.jt)' 


CO 


[6(co^j  + co) 


6{co^.j  - co)] 


(4.97a) 


2 

ne  X 


^ 1 + (<ot)2 


(4.97b) 


where  T f.  = N.,  n = N„/V,  and  co.  . = E.  ./Ti  (Q.E.D.) 

w 1 G G IJ  U 


ThG  CO  oo  Limit 

For  high  frequencies,  we  also  expect  to  recover  the  classical 
result  for  e.|(co)  because  the  quasi-continuum  limit  sets  in.  Thus, 
from  Equation  (4.44),  we  have 


e^(co) 


= 1 


2 

+ 

2 

CO 


1 Sire  V 47  r3 

^ h ' 


,r.j.|^  + 0(co'^)  (4.98) 


for  the  high  frequency  limit.  Therefore, 


lim  e.j  (co)  - 1 - — ^ 
or^  CO 


(4.99) 


as  expected. 


133 


The  Small  Particle  Problem 
Anomalous  FIR  Absorption 

The  FIR  absorption  coefficients  of  a = 20  A A1  and  Au  small 
particles  calculated  from  our  small  particle  model  are  shown  in  Figure 
4.7  for  three  different  ensembles.  Notice  that  our  calculation  clearly 
shows  the  effect  of  energy  gap  distributed  by  the  level  distribution 
function  on  the  FIR  absorption  coefficient.  In  other  words,  the 
"visibility"  of  the  quantum  size  effect  in  the  FIR  frequency  is 
enhanced  due  to  the  level  distribution.  Also,  the  effect  of 
magnetic  field  is  found  to  be  observable  as  a result  of  breaking  the 
time  reversal  symmetry  which  changes  the  ensemble  of  energy  levels 
from  orthogonal  to  unitary  ensemble. 

In  actual  experimental  situations,  there  exists  a particle 
size  distribution^®  as  an  inevitable  factor,  which  is  expected  to 
broaden  the  measured  absorption  coefficient.  The  effect  of  size 
distribution  on  the  FIR  absorption  coefficient  is  shown  in  Figure  4.8, 
along  with  the  absorption  coefficient  of  individual  ensemble  members. 
Also,  the  effect  of  the  size  distribution  on  the  FIR  absorption  for 
three  different  ensembles  is  shown  in  Figure  4.9,  while  the  detailed 
particle  size  dependence  of  the  FIR  absorption  coefficient  is  shown  in 
Figures  4.10,  4.11,  and  4.12  for  each  ensemble.  Notice  that  the 
frequency  dependence  of  the  FIR  absorption  coefficient  is  not  quadratic 
any  more  for  the  frequencies  away  from  the  gap  frequency  because  02(0)) 
defined  as  02(00)  = (1  - wei(a))}/4T:  becomes  much  smaller  than  o-jCw)  in 
magnitude  as  frequency  increases  above  the  gap  unlike  the  GE  case  and 


f-0.01  20A  A1  DorticlQS  In  vacuum  F“0. 01  20A  Au  particlQS  in  vocuum 

o.?r ' ' I I 
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Figure  4,7  FIR  absorption  coefficients  of  20  A size  A1  and  Au  particles  for  three  different 
ensembles 
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Figure  4.9  Effect  of  the  size  distribution  on  each  different  ensemble 
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Figure  4.10  Size  dependence  of  the  FIR  absorption  coefficient  for  Orthogonal  ensemble 
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Figure  4.12  Size  dependence  of  the  FIR  absorption  coefficient  for  Symplectic  ensemble 
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the  absorption  coefficient  is  sensitive  to  both  a-]  (u)  and  02(w)  as  can 
be  seen  from 


a 


e 


■V  fe 


3/2 

0 


2 


(4.100) 


and 


'v  feJ/2  Oi  a^  (0^  (4.101) 

where  (a^^)  is  the  classical  result  for  the  electric  (magnetic) 

dipole  obsorption  coefficient  discussed  in  Chapter  I. 

However,  the  quadratic  nature  in  frequency  dependence  in  a-|(a)) 

1 /3 

IS  expected  to  change  to  w dependence  above  the  gap  like  the  normal 

Anderson  insulator  above  the  crossover  frequency  (see  Appendix  E).  This 

effect  can  be  put  in  phenomenologically  by  introducing  the  semiclassical 
49 

approach  with  a proper  prescription  because  introducing  arbitrary 
broadening  does  not  conserve  the  local  electron  number. 

For  large  size  particles,  the  magnetic  dipole  contribution  can 
be  important  as  discussed  in  Chapter  I.  But,  Equation  (4.101)  totally 
depends  on  a classical  argument.  Therefore,  it  is  necessary  to  go  to 
the  multipole  calculation  beyond  the  dipole  approximation  in  order  to 
have  the  quantum  mechanical  justification.  However,  if  we  use  Equation 
(4.101)  to  describe  higher  order  effect  (5,  = 1 ) as  a good 

4 

approximation,  we  have  additional  co  contribution  to  the  absorption  as 
shown  in  Figures  4.13,  4.14,  and  4.15  and  this  additional  contribution 
changes  the  shape  of  the  absorption  coefficient  significantly  for  large 
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Size  dependence  of  the  FIR  absorption  coefficient  including  the  magnetic  dipole 
contribution  for  the  geometrical  standard  deviation  1.2 
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Figure  4.15  Size  dependence  of  the  FIR  absorption  coefficient  including  the  magnetic  dipole 
contribution  for  the  geometrical  standard  deviation  1.4 
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size  particles.  But,  the  magnitude  of  absorption  for  large  particles 
is  still  much  smaller  than  the  experimental  results  shown  in  Chapter 
III  and  the  frequency  dependence  is  not  correct.  This  discrepancy 
is  perhaps  due  to  the  incompleteness  of  the  theory.  Especially,  the 
approximation  given  by  Equation  (4.37)  requires  more  theoretical 
speculation  because  we  totally  excluded  the  frequency  dependence  on  the 
position  matrix  elements  in  the  approximation. 

Therefore,  we  will  restrict  the  theoretical  speculation  on  the 
FIR  problem  for  frequencies  much  smaller  than  the  cutoff  frequency, 
Vp/a.  Now,  for  small  frequencies,  we  can  examine  purely  particle  size 
independent  response  by  rescaling  Xg(w)  and  a-|(a))  given  by  Equations 
(4.59)  and  (4.68),  respectively,  as 


Therefore,  we  have  a scaled  FIR  absorption  coefficient,  which  is 
particle  size  independent  as 


Xe(“)  = Xg(w)/X 


0 

e 


(4.102) 


and 


(w)  = a-j  (w)/x° 


(4.103) 


(4.104) 


and  the  graphs  of  a(o))  for  the  three  different  ensembles  are  shown  in 
Figure  4.16. 
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Figure  4.16  Normalized  FIR  absorption  coefficient  to  three  different  ensembles 
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In  summary,  the  "anomalous"  FIR  absorption  is  truly  a quantum 
size  effect  resulting  from  the  broken  symmetry.  And  this  effect  can  be 
observed  experimentally  in  FIR  frequencies  for  small  particles  and 
microwave  frequencies  for  large  particles  (i.e.,  the  experimental 
test  purely  depends  on  the  size  quantization  gap).  However,  the 
extension  to  higher  frequencies  to  explain  the  experimental  results 
for  an  order  of  100  K particles  fails  mainly  because  of  the  incomplete 
approximation  in  the  matrix  element  calculation. 

Surface  Plasmon  Broadening 

As  shown  in  Table  1.1  in  Chapter  I,  a significant  broadening 

in  the  Maxwell -Garnett  (or  Mie)  resonance  peak  has  been  observed 

experimentally  which  does  not  count  on  either  classical  calculation  or 

Kawabata-Kubo  calculation.  In  addition,  a weak  temperature 

84 

dependence  of  the  resonance  strength  and  width  which  does  not  seem 
to  arise  obviously  in  classical  theory. 

A general  description  of  the  Maxwel 1 -Garnett  resonance  can  be 
obtained  from  Equation  (1.23)  as 

a(m)  = f ^ Im  Ye(“)  (4.105) 

Here  again,  a is  the  absorption  coefficient,  f is  the  volume  fraction 
of  small  particles  in  a medium  with  Cq,  and  is  the  electric 
polarizability  per  unit  volume  of  a particle  in  a medium,  e^. 

By  separating  the  dielectric  function  of  a small  particle, 
e|^(co)  into  real  and  imaginary  parts  as  e^(w)  = e-]  + ^^2* 
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Im  Ye(<»)  = 


:o  =2 


(e- 


2eo)‘ 


(4.106) 


from  Equation  (1.10).  Therefore,  the  absorption  coefficient  becomes 


a{o))  -f-y  eg 


(f 


2co)‘ 


2 

-2 


(4.107) 


Note  that  Equation  (1.27)  is  the  limiting  case  of  Equation  (4.107)  in 
the  limit  of  Cq  « e-j  for  the  electric  dipole  contribution  only. 

Now,  in  the  vicinity  of  the  resonance  frequency  w^p,  which 
satisfies  ^1^“^  expanded  as 


a-|  (w) 


Et  (co  ) 
1 ' sp^ 


de-j  (oj) 


do) 


(w 


CO  ) 

sp^ 


(4.108) 


sp 


because  e-j(co)  varies  slowly  at  high  frequencies.  Thus,  we  obtain 
a(co)  near  the  resonance  as 


^2<“sp> 


A‘(a.  - Ujp)^  + 


(4.109) 


A E 


de-|  (co) 
dco 


co=co 


sp 


where 
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From  Equation  (4,109),  we  can  clearly  see  that  the  line  width  is 
determined  by  £2^“)  ^ we,  therefore,  have  the  full  width 
at  half  maximum,  Aco,  as 


Aco  = 2|co^^2  " “spl 


2e2Kp^ 

A 


(4.110) 


Here  <^-^12  stands  for  the  frequency  at  half  maximum.  For  example.  Aw 
of  a Drude-like  small  particle  becomes 


<*“>Drude 


1 


Kp^>' 


(4.111) 


a 


as  found  elsewhere,^^’^®’^^ 

However,  the  quantum  mechanical  calculation  of  Aco  for  a 
48  90  49  91  92 

sphere  ’ or  a cube  ’ ’ with  the  smooth  surface  gives  Aco  close  to 

the  classical  result  within  less  than  20%  variation  except  for  the 

48  1 

result  of  Kawabata  and  Kubo,  which  gives  9 x , . Therefore, 

L.  C 1 d S S 1 Cd  1 

as  briefly  discussed  before,  the  role  of  the  surface  roughness  should 
be  taken  into  account  in  order  to  explain  the  broadening  observed  in 
this  work  and  in  the  work  of  others. 


CHAPTER  V 

SUMMARY  AND  CONCLUSION 


Since  the  anomalous  FIR  absorption  by  small  particles  was 
13 

discovered  by  Tanner,  many  efforts  have  been  made  to  explain  the 
peculiar  phenomenon  starting  from  how  to  take  the  average  of  the 
composite  system  (metal-insulator)  and  ending  with  clustering  theories. 
In  this  dissertation,  two  major  issues  have  been  raised  and  studied 
extensively:  One  is  the  clustering  theory  and  the  other  is  the 

possibility  of  nonmetallic  behavior  of  small  particles. 

For  the  first  issue,  we  have  demonstrated  the  experimental 
results  of  clustered  and  nonclustered  small  particle  composite  systems 
via  direct  FIR  and  MIR  absorption  measurements  and  the  TEM  studies 
and  we  have  verified  that  the  clustering  trfect  does  enhance  the  FIR 
absorption,  but  the  anomalous  FIR  absorption  observed  in  this 
experiment  is  truly  a small  particle  phenomenon!  In  addition, 
unexpected  new  observations  were  made  in  the  MIR  frequency  range, 
oxide  contribution  and  substantially  smaller  MIR  absorption 
coefficient.  Also,  the  grinding  effect  study  in  the  MIR  range  showed 
a clear  evolution  process  from  a clustered  to  a nonclustered  system, 
which  eventually  shows  the  linear  dependence  in  the  volume  fraction, 
f,  on  the  absorption  coefficient. 

A test  for  the  constant  relaxation  time  of  electrons  in  a 
small  particle  was  made  in  order  to  verify  the  second  issue,  nonmetallic 
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possibility.  Surprisingly,  what  we  observed  was  the  behavior  of  an 
insulator!  In  other  words,  the  intrinsically  poor  conductivity  of  a 
small  particle  is  responsible  for  the  puzzling  anomalous  FIR 
absorption  problem.  This  point  stimulated  the  theoretical  speculation 
on  the  fundamental  properties  of  a small  particle  and  we  found  the 
behavior  of  an  insulator  as  discussed  in  Chapter  IV.  However,  this 
theoretical  attempt  is  far  from  the  reality  to  explain  all  physical 
phenomena  of  small  particles,  which  we  have  observed  experimentally. 
Even  so,  we  believe  that  our  model  calculation  explains  the  correct 
behavior  of  a small  particle  near  the  gap  frequency  and  we  are  on  the 
right  track.  Hence,  further  extensive  theoretical  research  needs  to 
be  done  in  the  future  to  obtain  the  correct  description  of  a small 
particle  at  high  frequencies,  which  will  also  lead  to  the  answer  to 
the  surface  plasmon  broadening  problem. 

The  observation  of  the  red-shifted  interband  transition  and 
the  possible  "localized  band"  in  aluminum  small  particles,  which 
disappear  in  the  bulk  limit  as  shown  in  Appendix  D,  adds  up  to  a new 
feature  to  the  small  particle  problem. 

In  conclusion,  the  anomalous  FIR  absorption  is  arising  from 
the  intrinsic  properties  of  a small  particle,  which  cannot  be  explained 
by  the  classical  theory  and  remains,  therefore,  unsolved.  Also,  a 
small  particle  is  an  insulator  and  it  has  a different  energy  band 
structure  from  its  bulk  state. 
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APPENDIX  A 
SURFACE  ROUGHNESS 


As  described  earlier  in  Chapter  I,  the  roughness  of  the  surface 
of  a small  particle  to  atomic  scale  plays  a crucial  role  in  breaking 
degeneracies  of  eigenstates  and  alters  the  static  and  dynamic 
properties  of  a small  particle  significantly. 

If  we  model  a small  particle  as  a system  of  free  electron  gas 
confined  in  a sphere  and  treat  the  surface  as  a smooth  and  infinite 
potential  barrier,  then  we  can  get  all  necessary  information,  in 
principle,  by  solving  the  Schrbdinger  equation. 

H4>n(r)  = 


where  the  solution  4>p(>")  vanishes  at  the  boundary.  Therefore,  the 
corresponding  eigenvalues  and  eigenfunctions  are  found  as  follows: 
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(A, 2) 


and 


n£  a 


(A. 3) 


where  a is  the  radius  of  the  sphere,  is  the 

Bessel's  function,  k^^  is  the  zeros  of  j^(r),  and  c^^ 
normalization  constant. 


order  spherical 
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Now,  if  we  introduce  the  roughness  to  the  smooth  surface 


H(r')  'f^(r')  = E^'i'^(r') 


(A. 4) 


and  T^Cr')  = 0 at  the  boundary.  Here  r'  represents  another  coordinate 

for  a new  sphere  with  rough  surface.  If  the  degree  of  the  roughness 

is  not  serious,  we  can  treat  the  deviation  from  the  perfect  surface  as 
93 

a perturbation.  By  considering  one  parameter  mapping  r ' ->-r  = r(r' ,X) 
where  r is  an  analytic  function  of  X and  X is  the  mapping  parameter. 
Then,  r can  be  expanded  as 
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and  Jacobian  (J)  of  the  mapping  becomes 


9r 


9r' 


= 1 + X 


9J 


9X 


, 2 

f 2,1 

X 

9 J 

2 

. 2 

da 

+ O(x^) 


(A. 6) 


Thus,  the  Hamiltonian  H(r')  can  be  expressed  as  follows: 
H(r')  = H'(r,X) 


= H(r)  + X 


f9H'] 

9^H' 

[9Xj 

0 2 

9X^^ 

+ O(X^) 


(A. 7) 


with  H' (r,o)  = H(r). 

Now,  from  Equation  (A. 4),  a new  Hamiltonian  can  be  obtained 


such  as 
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H(r)  'F^(r)  = E^'F^(r) 


where  H(r)  = H(r')J^^^  and  ’l'^('^)  = 'l'^('^').  Next, 

Equations  (A. 6)  and  (A. 7),  we  have 


and 
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H(r)  = H(r)  + XH^  + T ^2  ^ 
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Therefore,  by  expanding  in  powers  of  X, 


= <()^(r)  + Y + o(X^) 


and  defining  a similar  expansion  for  E^, 


E„  = e„  + AE„ 
n n n 


(1), 


+ O(x^) 


(A. 8) 

from 


(A. 9a) 
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(A. 10) 


(A. 11) 
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we  can  apply  a standard  perturbation  theory  to  get  the  correction 
terms. 

For  a simple  case  r'  = r(l  + X),  an  easy  justification  of  the 

perturbation  treatment  can  be  made  by  solving  Equation  (A. 4)  with  the 

condition  H'  [a(l  + X)1  = 0.  By  considering  transformations, 
n m 


6 = 0 


(A. 12) 


<p  = (p 


then,  the  Jacobian  for  the  transformation  becomes 


J = 


1 + X 


= 1 - X + X^  + . . . 


(A. 13) 


Hence, 


H = J‘^^^  Hf(l  + X)  r,  0,  (p\  J^/^ 


= H(r,  6,  (|))  - 2X  H(r,  6,  <j))  + 3X^  H(r,  6,  <})) 


(A. 14) 


and 
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(A. 15) 
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Thus,  we  obtain 


XE 


= -2Xe  — = 3X^e 


(A. 16) 


where  is  given  in  Equation  (A. 2).  And 
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due  to  the  orthogonality  of  the  unperturbed  eigenfunction  (J>^  given  in 
Equation  (A. 3). 

Now,  let's  consider  the  exact  solution.  The  eigenfunction 
TCr ' ) can  be  written  as 


¥(r' ) 
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with  eigenvalue 
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From  Equation  (A. 8), 


'F^(r')  = ^^(r) 


(A. 21) 


therefore. 
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which  agrees  with  Equation  (A. 19)  to  the  second  order  with 


I 2^8^ 


(A. 23) 
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Tavel  carried  out  a detailed  perturbation  calculation  for  a mapping 


r'  = r[l  + Xf(e,  (fi)] 


(A. 24) 


e'  = 0 


<f>'  = 4> 


and  by  expressing  the  surface’ perturbation  Xf(0,  <}>)  in  the  most  general 
form 


Xf(e,  (j))  = X 


1 I 

£=1  m=-S, 


+ Y„  ) + 
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(A. 25) 


where  Y^  are  the  spherical  harmonics  and  each  of  the  coefficients  a^ 

and  b.„  is  chosen  from  uncorrelated  Gaussian  distributions  with  zero 
£m 
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mean  and  unit  standard  deviation,  following  important  results  are 
found  in  this  model,  after  a tremendous  amount  of  calculation: 

1.  There  is  an  upper  limit  to  the  value  of  Z which  can 

break  the  degeneracy  of  any  given  level.  In  this  model, 

SL  =13.  Therefore,  each  surface  in  the  ensemble 
max 

gives  rise  to  a 27  x 27  matrix  and  only  spherical 

i U 

harmonics  up  to  the  26^^  order  in  the  surface 
perturbation  contribute  to  the  level  splitting  from  the 
Wigner-Eckart  theorem. 

2.  The  degree  of  roughness  can  be  varied  by  multiplying 
the  harmonic  coefficients  by  inverse  powers  of  Z to 
damp  out  the  contributions  of  higher  order  spherical 
harmonics.  For  the  maximally  rough  case  (Z  ^ damping 
with  X = 0.6),  a distribution  of  level  spacings  which 
agrees  with  the  Wigner  distribution  is  obtained,  thus 
supporting  the  applicability  of  the  random  matrix 
theory. 

3.  There  is  a significant  correlation  between  the 
eigenvalues,  regardless  of  their  particular  symmetry, 
even  though  the  matrix  element  themselves  are  chosen  to 
be  uncorrelated.  This  correlation  reflects  the 
consequence  of  the  level  repulsion  predicted  by  Wigner. 

4.  As  the  degree  of  the  roughness  decreases,  a significant 
change  occurs  in  the  level  distribution  which  is  rather 
close  to  the  level  attraction  case. 
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Although  this  model  calculation  is  an  idealized  case  and  does 
not  describe  the  realistic  small  particle  surface,  the  idea  of 
correlation  arising  from  the  random  irregularity  of  the  surface 
would  be  feasible  and  there  is  no  doubt  about  the  necessity  of  new 
classes  of  ensemble  which  can  describe  the  level  correlation  and  its 
distribution  in  a small  particle. 


APPENDIX  B 
LEVEL  STATISTICS 


4 

Kubo  suggested  that  the  effect  of  random  surface 
irregularities  of  a small  particle  randomize  the  distribution  of  energy 
levels  and  the  author  took  the  adjacent  energy  gap  distribution  to  be 
a Poisson  distribution 


P(e) 


e-s/AE 

AE 


(B.l) 


where  e is  the  energy  level  spacing  and  AE  is  the  average  spacing. 

Gor’kov  and  El iashberg^  discussed  the  effect  of  level  repulsion 

between  two  crossing  eigenvalues  by  employing  the  random  matrix  theory 
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introduced  first  by  Wigner  and  developed  further  by  Dyson,  which  had 

proven  to  be  quite  realistic  in  describing  excited  nuclear  states. 

Because  the  eigenvalues  of  the  nuclear  Hamiltonian  are 

essentially  unknown  due  to  the  extensive  energy  level  configuration 

mixing  and  there  exists  a well-established  connection  between  the 

spacing  distribution  of  adjacent  eigenvalues  of  a matrix  with  random 

chosen  elements  and  the  level  spacing  distribution  of  excited  nuclear 

states,  it  is  a logical  conclusion  to  make  the  same  analogy  to  a small 

particle  whose  Hamiltonian  describing  electrons  in  the  particle  is  not 

known  due  to  some  complex  interactions  between  electrons  and  the 

impurities  and  irregularities  of  the  surface. 
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Wigner  first  made  a statistical  hypothesis  that  all  matrix 
elements  of  a system  having  N energy  levels  are  statistically 
independent  and  satisfies  the  equal  a priori  condition.  However,  one 
cannot  define  a uniform  probability  distribution  on  an  infinite  range 
and  this  problem  can  be  avoided  by  introducing  a cutoff  into  the 
theory  which  could  correspond  to  excluding  the  strong  interactions  in 
the  nuclear  theory.  The  best  description  of  such  a system  is  the 
Gaussian  distribution 


N 

H 

P(H)  = A exp 

- 1 

mn 

2 

m,n 

where  A is  the  cutoff  excluding  the  divergence. 

Later,  Dyson  conjectured  that  it  is  impossible  to  define  an 
ensemble  in  terms  of  the  H in  which  all  interactions  are  equally 
probable  and  the  requirement  of  the  cutoff  was  imposed  without  clear 
physical  motivation.  Also,  the  author  pointed  out  that  it  is  not 
reasonable  to  picture  the  as  resulting  from  some  random  process  of 
a conventional  statistical  ensemble. 

In  the  phenomenological  Dyson  theory, the  Hamiltonian  matrix 
element  description  was  rejected  and  a unitary  N x N matrix  S which  is 
related  to  the  Hamiltonian  of  the  system  in  a complicated  way  was 
considered.  According  to  Dyson's  hypothesis,  the  correlation  of  n 
consecutive  energy  levels  of  a system  (n  « N)  is  statistically 
equivalent  to  that  in  the  ensemble  of  n consecutive  angles  6^  on  the 
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unit  circle.  The  connection  between  S and  H is  left  to  be  vague  to 
reality,  but  a definite  relation,  for  example,  may  be  considered  such 
as 


and  S can  be  diagonalized  in  terms  of  a unitary  matrix  U as  follows: 


Now,  the  unitary  matrix  S can  have  a definite  mathematical 
symmetry  which  depends  on  the  existing  symmetries  in  the  real  system. 

If  the  Hamiltonians  have  timer  reversal  and  spin  rotation 
invariant  symmetry,  then  all  wave  function  can  be  made  to  be  real. 
Thus,  all  the  matrix  elements  of  the  Hamiltonian  are  real  in  this  call 
and  an  ensemble  of  such  Hamiltonian  is  invariant  under  orthozonal 
transformation  (i.e.,  Hamiltonian  of  a system  with  even  electron  spins 
and  zero  external  magnetic  field).  The  matrix  S has  the  same  symmetry 
as  H which  is  obvious  from  Equation  (B.3).  Therefore,  now  Equation 
(B.4)  can  be  made  to  be  real  and  U takes  the  form 


S = exp  [iH],  S = 


(B.3) 


S°  = U'^SU 


(B.4) 


ie. 

with  (S°)..  = 6..  e ^ and  U U = 1. 

I O I «J 


(B.5) 


where  U is  a real  matrix. 


When  the  time  reversal  invariant  symmetry  is  broken,  the  matrix 
elements  of  H can  no  longer  be  real.  Then  the  unitary  matrix  S does  not 


163 


have  any  additional  symmetry.  Therefore,  U is  now  an  arbitrary 
unitary  matrix  and  this  symmetry  is  suitable  for  the  system  in  the 
presence  of  a magnetic  field.  The  third  possible  class  of  symmetry 
is  the  simplectic  case  that  a system  has  the  time  reversal  symmetry, 
but  has  no  rotational  symmetry  (i.e.,  a system  with  odd  number  electron 
spins  in  a zero  magnet  field  or  one-electron  Hamiltonian  with  strong 
spin  orbit  coupling  in  a zero  magnetic  field).  Thus,  all  levels  in 
the  system  have  Kramer's  doublets  due  to  spin-orbit  interactions 
(i.e.,  the  spin  and  the  orbital  states  are  inseparably  coupled 
together),  and  each  eigenvalue  of  S has  to  appear  twice. 

Now  the  concept  of  equal  probabilities  for  all  matrix  elements 

of  S in  calculating  the  average  implies  an  equal  contribution  of  each 

elementary  volume  AS  on  the  hypersurface  {S}.  In  terms  of  the 

2 

coordinates  6^-  and  U,  the  elementary  area  (dL)  becomes 

(dL)^  = Tr(dSdS'^)  (B.6) 


Now,  from 


dS  = UV6U  + 6U‘^S°U  + U'^(ide.)U  (B.7) 

we  obtain, 

(dL)^  = 2Tr  (S°6U6UV‘^  - 6US°'^6UV)  + \ (de.)^  (B.8) 

where  the  property  6UU^  = -U6U^  have  been  used.  Equation  (B.8)  can  also 
be  written  as 
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(dL)2  = I (6U,j6U;j) 
1 »J 


2 sin  (0^. 


+ I (de.)' 


(B.9) 


by  using 


^ J 


Thus, 
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The  volume  element  AS  can  be  calculated  by  considering  the  off  diagonal 
elements  such  as 


AS 


orth 
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i<j 


n de^AU 


(B.ll) 


where  AU  is  the  volume  element  of  vector  <5U...  Therefore,  the 

' V 

probability  of  finding  a system  of  orthozonal  group  belongs  to  the 
volume  of  AS  is 


p(AS)  = j AS 


(B.12) 


where  V = 


AS,  is  the  total  volume.  For  unitary  the  ensemble  case. 


since  all  the  elements  of  the  Hamiltonian  are  not  real,  there  are  three 


independent  variables  for  two  adjacent  level  Hamiltonians  by  choosing 
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the  off-diagonal  term  a-|2  to  be  a-j2  + Consequently,  we  have  to 

consider  a transformation 

U^GU  (B.13) 

i^i 

if  G is  the  diagonal  unitary  matrix  with  (G)^.  = e where  are 
an  arbitrary  angle  on  a unit  circle.  Therefore,  the  total  volume 
element  in  this  group  will  be 
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unit 


= AS 


orth 


AG 
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thus. 
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n de^.AU  (B.15) 
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Now,  for  the  symplectic  ensemble  case,  we  have  a Hamiltonian  consisting 
of  quarternoions . In  other  words,  we  have  two  more  random  variables  in 
the  eigenvalue  of  the  Hamiltonian;  therefore,  two  more  volume  elements 
have  to  be  considered  because  the  unitary  matrix  U is  also  symplectic. 


AS 


unit 
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with  AG  = n dx^,  a = 1,  2,  3. 
a,j  ^ 
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A constant  density  of  states  can  be  calculated  from  the 

probability  distribution  because  it  depends  only  on  the  differences  of 

0.  and  0.. 

• 0 

The  next  interesting  quantity  is  the  two-level  correlation 
function,  which  determines  the  probability  that  a level  falls  in  an 
energy  interval  e from  another  level.  Dyson^  calculated  the  following 
correlation  functions  for  three  different  ensembles: 

„ort  . lljA  _ 
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for  the  orthogonal  ensemble. 


(x) 


1 - 


. 2 
sin  X 


(B.18) 


for  the  unitary  ensemble,  and 

RSympl  (x)  = 1 - + 

for  the  symplectic  ensemble.  Here  x = ^ , AE  is  the  average  level 
spacing.  The  detailed  behavior  of  R(x)  is  shown  in  Figure  B.l. 

Notice  that  the  two-level  correlation  function  has  essentially  the 
same  behavior  as  the  pair  correlation  function  of  identical  hard 
spheres  with  the  radius  AE.  However,  these  correlation  functions  have 
completely  different  limiting  cases^  as 
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Figure  B.l  Two  level  correlation  function  for  three  different 
ensembles 
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APPENDIX  C 
INTERBAND  TRANSITION 

The  optical  properties  of  simple  polyvalent  metals  are  known 
to  show  an  interband  transition.  Because  these  nearly  free  electron 
metals  have  the  inverse  square  root  singularities  in  the  joint  density 
of  states,  which  has  direct  relation  with  the  otpical  conductivity, 
these  metals  have  the  optical  absorption  associated  with  the  energy 
gap  arising  from  the  singularity. 

Theoretically,  the  concept  of  pseudopotential,  which  is  on 
the  basis  of  the  idea  of  the  screening  of  the  ion  by  electrons,  allows 
us  to  construct  a model  of  nearly  free  electrons,  which  now  feels  only 
the  weak  pseudopotential.  However,  unlike  the  free  electron  case,  the 
nearly  free  electrons  suffer  the  coherent  scattering  by  the  Bragg  plane 
even  in  the  weak  pseudopotential  which  greatly  modifies  the  optical 
properties  of  metals.  In  other  words,  the  existence  of  a periodic 
pseudopotential  introduces  an  energy  gap  at  Bragg  planes.  Thus,  when 
the  wave  vector  t.  of  an  electron  varies  from  one  branch  to  another 
continuously,  the  energy  changes  discontinuously  at  ^/2  by  2jU^j  when 
the  periodic  pseudopotential  is  expanded  in  a Fourier  series  such  as 

U(f)  = I Uj  e’®-"  (C.l) 
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where  ^ is  the  reciprocal  lattice  vector.  Now  can  be  determined  by 
the  deHaas-van  Alphen  effect  analysis  or  by  the  optical  measurements. 

A pseudopotential  calculation  of  the  aluminum  energy  band  done 
by  Ashcroft  found  that  a pair  of  bands  may  be  substantially  parallel 
when  they  are  plotted  in  a plane  parallel  to  the  zone  face.  By  using 
this  fact,  there  are  two  different  interband  transitions:  one  is  the 

parallel  band  transition  and  the  other  is  the  normal  interband 
transition. 

The  electron  wave  function  away  from  the  Bragg  plane  can  be 

regarded  as  a plane  wave  of  a free  particle,  and  the  wave  function  near 

the  Bragg  plane  can  be  written  as  a superposition  of  two  plane  waves. 

98 

Ashcroft  and  Sturm  carried  out  the  interband  conductivity  in  a local 
pseudopotential  approximation,  and  found  the  following  results  by 
separating  the  lower  (-)  and  upper  (+)  bands: 
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with 


(t  - 2)2  - t2 
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and  the  energy  in  the  upper  bands  (+)  and  lower  bands  (-)  as 


IyI  ± (1  + 


(C.5) 


where  k has  been  resolved  into  its  components  parallel  (k||)> 
perpendicular  (tp  to  and  y is  a function  of  k||  only  by  definition. 

Now,  as  shown  in  Figure  C.l(a),  the  transition  for  the  lower 
band  to  the  upper  band  has  a constant  interband  energy  difference, 
which  is  called  "parallel  band"  absorption.  Also  the  normal  interband 
absorption  is  shown  in  Figure  C.l(b).  Clearly,  these  two  different 
contributions  are  categorized  as  follows.  For  parallel  band 
absorption,  the  condition  written  below  should  be  satisfied  at  the 
Fermi  energy 


There  is  a high  frequency  cutoff  in  the  parallel  band  transition 
which  comes  from  the  eigenvalue  equation  for  kp  at  E = Ep 


(C.6) 


Therefore,  we  have 


(C.7) 


(C.8) 
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Figure  C.l  Aluminum  energy  bands  (from  Reference  98) 

(a)  Plotted  for  k near  the  L point  zone  face  where 
parallel  band  character  is  exhibited  and  (b)  normal 
interband  absorption 
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or 
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(C.9) 


to  the  first  order  in  (U^/Ep).  And  the  gap  energy  becomes 


- E_  = 2 


U^l  (1  + 


(C.IO) 


from  Equation  (C.5)  or  (C.7). 

When  the  upper  band  is  above  the  Fermi  energy  and  the  lower  band 

is  below 
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and 
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then  we  are  in  the  region  of  the  normal  interband  transition  as  shown  in 
Figure  C.l(b).  Thus,  we  have 


{k|}^  - {kp^  = 


2 lU: 


[y  - (1  + 


^2 

(C.12) 


Equation  (C.12)  can  be  rewritten  as 


{kl}^  - {k|}^  - (710)  + 7io)(,)  (Tio),  - Tito) 

^ ^ 4e4  u I 


(C.13) 
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where 


hco  = 2|U^|  (1  + 


and 


flco,  = 2(Eg£p  . u2)’/2  . Ef 

Therefore,  there  is  an  upper  cutoff  frequency  for  the  normal 
interband  transition  with  tuo-i  = Tioa 
place  between  TiaiQ  and  tuo^ . 

For  aluminum,  there  are  two  parallel  transitions  at  U-iii  and 
^200’  ^^owever,  the  pseudopotential  U-iii  gives  an  absorption  at 
0.44  eV  (3,550  cm~^)  which  is  weaker  than  the  absorption  dut  to  U2QQ. 
U200  ''5  responsible  for  the  absorption  at  1.52  eV  (=  12,260  cm  ^), 
which  is  strong  enough  to  be  seen  in  ordinary  optical  relection 
measurement. 

One  remark  on  this  type  of  calculation  is  that  a broadening 

effect  needs  to  be  considered  for  a real  system  as  done  by  Ashcroft 
98  99 

and  Sturm.  But  Benbow  and  Lynch  found  a good  agreement  with  the 
Ashcroft  and  Sturm  theory  after  introducing  separate  relaxation  time 
for  0.44  eV,  and  Drude  contributions  to  the  optical  conductivity. 


g + 2E^  and  the  transition  takes 


APPENDIX  D 

REVOVERY  OF  THE  BULK  NATURE 

As  we  have  seen  in  Chapter  III,  another  peculiar  behavior  of 
the  small  particle  is  the  interband  transition  in  small  aluminum 
particles.  If  this  peculiar  behavior  is  a true  small  particle 
phenomenon,  then  these  characteristics  will  be  lost  when  these  small 
particles  recover  their  bulk  characteristics.  One  idea  is  to  make  a 
pure  aluminum  pellet  by  compressing  the  aluminum  small  particles  and 
then  study  the  reflectivity  of  the  aluminum  pellet. 

The  measured  DC  conductivity  of  the  aluminum  pellet  was  in  the 
range  of  100-200  cm”\  which  is  a rather  low  value  to  judge  the 
complete  recovery  of  the  bulk  characteristics.  However,  since  the 
aluminum  pellets  are  slightly  over  the  percolation  threshold,  we  can, 
therefore,  expect  that  the  single  small  particle  characteristics  will 
disappear  to  some  extent  since  particles  are  touching  one  another  in 
this  regime. 

In  order  to  see  the  interband  transition  clearly,  a silver 
mirror  was  used  from  FIR  to  20,000  cm  ^ frequency  range.  Beyond 
20,000  cm”\  an  AL  mirror  was  used  and  normalized  to  an  Ag  mirror 
since  Ag  has  the  screened  plasma  frequency  at  3.9  eV  ('^31,500  cm  ^). 

The  reflection  measurements  for  two  different  aluminum  pellets 
are  shown  in  Figure  D.l;  one  was  made  with  <r>  = 103  ^ A1  small 

O 

particles  and  the  other  was  made  with  <r>  = 280  A particles.  There  are 
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Figure  D.l  Reflection  data  for  two  compressed  aluminum  small  particle  sample 
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two  noticeable  changes  from  Figure  B.17  in  the  plot:  One  is  the 

recovery  of  the  bulk  interband  transition  peak  and  the  other  is  the 
trace  of  the  broad  absorption  around  20,000  cm  \ which  was 
interpreted  as  a result  of  the  surface  localization  which  still  exists 
in  <r>  = 103  A particle  data,  but  not  in  <r>  = 280  A data.  Also,  the 
result  of  the  Kramers-Kronig  analysis  is  shown  in  Figure  0.2.  Note 
tha.t  the  overall  ac  conductivity  is  quite  similar  to  that  of  small 
particles  obtained  in  Chapter  III. 


Optical  conductivity  of  compressed  small  particle  pellet 
A1  14  <r>=103A 

A1  17  <r>=280A 
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Figure  D.2  Optical  conductivity  of  compressed  aluminum  small  particle  pellets 


APPENDIX  E 

AN  UNSUCCESSFUL  MODEL  CALCULATION 

Ordered-Disordered  Metal  Model 
As  discussed  in  Chapter  II,  small  particles  have  an  oxide 
coating  on  their  surface.  Since  the  electron  diffraction  study  of 
small  aluminum  particles  did  not  show  any  direct  evidence  of 
wel 1-crystal ized  aluminum  oxide  and  the  experimental  study  in  MIR 
frequency  done  in  Chapter  III  suggests  that  the  degree  of  amorphism 
in  the  oxide  coating  of  a small  aluminum  particle  is  even  higher  than 
the  result  of  Reference  68  (remember  that  the  amorphism  tends  to  push 
the  longitudinal  modes  of  aluminum  oxide  towards  the  higher 
frequencies),  it  has  been  assumed  that  the  oxide  was  formed  in  an 
amorphous  state  or  formed  in  a nonstoichiometrical  way.  However,  in 
actual  calculation,  we  have  assumed  that  the  amorphous  oxide  is 
loss-free  or  negligible,  if  any.  This  assumption  requires  somehow 
further  speculation  because  as  we  can  see  from  Equation  (4.100),  the 
FIR  absorption  can  be  enhanced  if  a-|  « *^D^Lide‘  Sievers^*^*^ 

observed  the  similar  anomalous  FIR  absorption  by  fully  oxidized 
particles  (evaporated  in  O2  atmosphere).  Additionally,  the 
experimental  evidence  on  the  insensitivity  of  the  absorption 
coefficient  to  the  particle  size  variation  suggests  that  the  observed 
anomalous  FIR  absorption  may  not  be  due  to  the  quantum  size  effect 
of  pure  electron  gas  confined  in  a restricted  volume,  rather  may  be  due 
to  the  amorphous  oxide  coating  with  some  finite  but  poor  dc 
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(a) 


Figure  E.l  A small  particle  model  incorporated  with  Anderson  insulator 
(a)  An  illustration  of  a small  particle  and  (b)  general 
behavior  of  Anderson  system  for  various  dimensions  (from 
Reference  101) 
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54 

conductivity  as  examined  by  Sen  and  Tanner  or  with  some  combination 
of  poor  ac  o-j  and  as  mentioned  above  if  the  amorphose  oxide  is  an 
insulator. 

Judging  from  this  evidence,  we  can  construct  a small  particle 
model  which  consists  of  two  parts  as  shown  in  Figure  E.l(a):  One  part 

is  the  core  part  with  crystalline  (ordered)  metal;  the  other  is  the 
outer  shell  formed  with  the  disordered  metal  (metal  with  lots  of 
impurities  and  lattice  defects),  which  we  view  within  the  Anderson 
localization  model  as  a disordered  electron  system, 

Anderson  Localization 

The  disorder  driven  metal-insulator  transition  which  is  known 
as  the  Anderson  transition  is  well  understood  and  generally  believed 
to  occur  in  a 3-D  system  but  not  in  2-D  and  1-D  systems.  By 
considering  a tight-binding  one-electron  Hamiltonian, 


I 


^ + 
E.c^Ci 


I Jij 


+ 

C.c. 

• VJ 


(E.l) 


where  E^.  is  the  energy  level  of  an  electron  at  site  i,  J..  is  the 
hopping  integral  between  site  i and  j,  and  c|  (c^. ) is  the  creation 
(annihil ization)  operator  of  an  electron  at  site  i.  Now,  if  we  carry 
out  a standard  tight  binding  expansion  of  the  eigenfunction'!'  with 
eigenvalue  E in  terms  of  atomic  orbitals  such  as 

'F  = I (E.2) 

i 
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where  is  the  atomic  orbital  centered  at  site  i.  Therefore,  we 
find  an  equation  for  the  expansion  coefficient  as 


Now,  in  the  case  of  the  nearest  neighbor  hopping  only,  we  have  an 
expression  for  the  coefficient  a^.  such  as 

E a.  = E.  a.  + I JQl^^i-l  ^i+1^  ^ 

where  z is  the  number  of  the  nearest  neighbors. 

Next,  we  want  to  set  up  a criterion  for  the  metal-insulator 
transition  in  terms  of  the  random  impurity  potential  height  w and 
the  band  width  in  the  absence  of  disorder,  U.  When  w = 0,  we  have  a 
completely  extended  regime  (i.e.,  E^  =0),  thus  we  find 

^ ^m  " I '^01  ^Vl  Vl^  ^ 

By  assuming  a plane  wave  solution  for  the  coefficient,  a^  as 


we  obtain  a solution 


E a, 


(E.3) 


E = z cosa 


(E.7) 


and  the  band  width  becomes 


U - 2z 


(E.8) 
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Next,  by  starting  from  U = 0 localized  limit  and  treating  the 
hopping  term  in  Equation  (E.l)  as  a perturbation,  we  can  determine 
the  localization  criterion  by  simply  examining  the  convergence 
condition  for  the  perturbation  series  expansion.  Assuming  a uniform 
energy  distribution  over  the  nearest  neighbor  sites  (i.e.,  w = 2z  |E^ 

- E . I ) , we  have 

<i  I H-|  I j>  2z  Jqi  y 

^ ' ' = — = - < 1,  the  condition  of  convergence 

E.  - E.  w w 

* J 

(E.9) 


from  the  amplitude  of  mixing  of  the  unperturbed  states  i with  its 
neighboring  states  j to  the  first  order.  In  other  words,  if  we  can 
treat  the  hopping  term  as  a perturbation,  electrons  are  localized. 
Otherwise,  electrons  are  in  the  extended  state.  Therefore,  we  have 
the  condition  U/W  < 1 for  the  Anderson  localization. 

Another  important  aspect  of  the  Anderson  localization  is  the 
scaling  behavior  discussed  by  Thouless^^  and  Abrahams  et  al.^^^  They 
considered  a block  of  solid  containing  many  sites  with  volume  in 
d-dimension  and  by  comparing  two  energy  scales  (a  similar  way  as 
done  in  Chapter  I),  the  dimensionless,  scale-dependent  conductance 
g(L)  = 6E(L)/AE(L)  was  obtained  as 

g(L)  = 

e 

2 

where  c?q  = e Dv(Ep)  and  D is  the  diffusion  constant.  The  key  point 
of  the  formulation  done  by  Thouless  is  that  once  we  know  a fundamental 
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conductance  gg  for  a system  with  a cube  of  side  Lg,  the  conductance 
g(bLg)  is  completely  determined  by  the  fundamental  conductance  gg  and 
the  scale  factor  b.  If  we  treat  b as  a continuous  transformation  such 
as 

g[(l  + e)L]  = f(e,  g(L)],  c « 1 (E.lla) 


= g(L) 


d in  gar 

^ d in  I 


(E.llb) 


then  we  obtain  the  scaling  behavior  which  can  be  specified  by  a 
differential  scaling  function  3(g): 


3(g)  = 


d £n  g(L) 
d in  L 


(E.12) 


From  Equations  (4.121)  and  (4.123),  the  limiting  cases  can  be  found  as 
1 im  3(g)  = d - 2 (E.13) 


and 


lim  3(g)  = ^.n  g (E.14) 
g-i-0 

The  detailed  behavior  of  function  3(g)  is  shown  in  Figure  E.l(b).  Note 
that  there  exists  an  unstable  fixed  point  g^  which  satisfies  3(g^)  = 0 
only  for  d = 3.  The  existence  of  g^  signals  the  Anderson  transition. 
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Conductivity  near  the  Anderson  Transition 

1 02 

Vollhardt  and  Wolfle  obtained  a scaling  function 


B[g(L)]  = d X-n  g/d  Jin  L which  is  consistent  with  the  results  given  in 
the  previous  section  by  using  the  self-consistent  diagrammatic  theory. 
Also,  they  found  the  generalized  diffusion  coefficient  D(co)  of  a 
disordered  system  within  the  same  context. 

The  frequency  dependent  conductivity  in  the  disordered  electron 
system  can  be  found  by  calculating  the  diffusion  coefficient  D(co)  for 
2 < d < 4 given  by 


disorder  parameter.  Here  the  cutoff  is  given  by  kp  = (irAkp). 
Rearranging  Equation  (E.15)  gives 


(E.15) 


1 » 

where  the  bare  diffusion  constant  Dq  = (frdmA)"  with  X = ^ , the 


+ Xd  — ^ 


(d-2)/2  rL(o))  d-3 

j A 


dx  

•'0  1 + X 


(E.16) 


where  L(6j)  - [D(a))/(-io)TDg)]^'*^^.  Thus,  we  have  the  critical  coupling 

constant  X from 
c 


X 


X 


c 


(E.17a) 
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Here 


X 


c 


d - 2 

fk  ^ 

z-a 

fd  - 2 2-dl 

Q. 

l^rj 

II 

Q_ 

l/(d-l) 


(E.17b) 


For  COT  « 1 - X/X^,  Equation  (E.16)  can  be  solveci  by  taking 
the  L(co)-x”  limit  in  the  conducting  regime  (X  < X^)  and  in  the  insulating 
regime  (X^  < X)  to  obtain 


g(to)  ^ 
^0 


(d-2)/2 


o{co)/aQ 


(E.18) 


where  = r(d/2)r(2-d/2) . 
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For  d = 3,  the  cubic  equation  can  be  solved  to  find 


oM . 

in  the  conducting  regime  (X 
the  solutions  are 

oM  = ( )V3 

Oq  C 


[1  + -^  (-i  /co^)^^^],  CO  « (0^ 


' ' 
IT 

2/3 

T-icol 

’/3  ^ 2l 

2 

- 

^ J 

i 

< X^).  In  the  insulating  regime 


f \ 

2 

f2 

2 • 1 

7T 

IT 

CO 

ICO 

2, 
V.  J 

2 

“c 

CO  « CO, 


* ^ 
TT 

2/3 

■ 

ico 

1/3  2^ 

2 

V.  J 

- 

" 3 

c 

4 

CO  » CO, 


(E.19a) 

(E.19b) 
(X  > X^,), 

(E.20a) 

(E.20b) 
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Here  the  crossover  frequency  is  defined  as 


(E.21) 


and  oj  also  can  be  written  as 
c 


(E.22) 


in  terms  of  the  localization  length 


FIR  Absorption  Coefficient 


Now,  the  FIR  absorption  coefficient  of  our  new  model  is  ready 
to  be  calculated  by  using  the  formula  given  by  Equation  (3.6).  Since 
we  are  dealing  with  fairly  large  particles  ('^^100  K)  in  the  experiment, 
we  describe  the  core  part  of  a small  particle  with  Drude  dielectric 
function  as  a good  approximation.  In  fact,  the  anomalous  FIR  absorption 
does  not  come  from  the  core  part  in  this  model.  Hence,  it  is  almost 
independent  of  the  choice  of  for  the  core  part  of  a small 

particle. 

For  the  oxide  coating,  a-|(a3)  and  02(w)  calculated  in  Equation 
(E.20)  is  used.  As  a result,  a remarkably  good  agreement  between  this 
model  and  the  FIR  experimental  result  was  obtained  as  shown  in  Figure 
E.2.  Notice  that  the  contribution  of  the  magnetic  dipole  term  in 
addition  to  the  electric  diple  term  shows  a nice  fit  for  large  size 
particles. 
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Figure  E.2  A fit  obtained  from  the  model  calculation 
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However,  this  model  calculation  fails  to  explain  the 

superconducting  state  behavior  of  small  particles  as  observed  by  Carr 
104 

et  al . unless  the  Anderson  insulator  undergoes  a superconducting 
transi tion. 
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